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This paper is dedicated to the memory of our friend and colleague John L. Norton, who 
wrote the original versions of the computer programs that we use to calculate the single-particle 
energies and resulting shell and pairing corrections for a deformed folded- Yukawa single-particle 
potential. 
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Abstract: We tabulate the atomic mass excesses and nuclear ground-state deformations of 8979 nuclei 
ranging from 16 to A = 339. The calculations are based on the finite-range droplet macroscopic 
model and the folded- Yukawa single-particle microscopic model. Relative to our 1981 mass table 
the current results are obtained with an improved macroscopic model, an improved pairing model 
with a new form for the effective-interaction pairing gap, and minimization of the ground-state 
energy with respect to additional shape degrees of freedom. The values of only 9 constants are de- 
termined directly from a least-squares adjustment to the ground-state masses of 1654 nuclei ranging 
from 16 to 263 106 and to 28 fission-barrier heights. The error of the mass model is 0.669 MeV for 
the entire region of nuclei considered, but is only 0.448 MeV for the region above N — 65. 



1 Introduction 

We presented our first macroscopic- microscopic global nuclear mass calculation 12 years ago 1,2 ). 
This calculation, which was based on a finite-range liquid-drop model for the macroscopic energy 
and a folded- Yukawa single-particle potential for the microscopic corrections, was somewhat 
limited in scope. With only 4023 nuclei included, it did not extend to the proton or neutron 
drip lines or to the region of superheavy nuclei. Also, the quantities tabulated were limited to 
ground-state masses, Q2 and Q4 moments, and microscopic corrections. 

Our next publication of calculated nuclear masses occurred five years ago 3 ' 4 ). In these 
calculations new pairing models had been incorporated and two different macroscopic models 
were investigated, namely the finite-range liquid-drop model (FRLDM) 3 ) and the finite-range 
droplet model (FRDM) 4 ). These abbreviations are also used to designate the full macroscopic- 
microscopic nuclear structure models based on the respective macroscopic models. The former is 
the macroscopic model used in the 1981 1 > 2 ) calculations and the latter is an improved version 5 ) of 
the droplet model 6 " 8 ). Because there were several unresolved issues in the 1988 calculations 3 ' 4 ) 
these tables should be regarded as interim progress reports. 

We have now resolved these issues, which were related to the pairing calculations 9 ), to the 
effect of higher- multipole distortions on the ground-state mass 10 ), and to some details of the 
shell-correction and zero-point-energy calculations 11 ). The resolution of these issues has resulted 
in the present mass table. We first briefly review some important results obtained in the original 
1981 calculation and enumerate the additional features of our new calculations. 

Subsequent comparisons of predictions of our original 1981 model 1 ' 2 ) with nuclear masses 
measured after the calculations were published showed that the model would reliably predict 
masses of nuclei that were not included in the determination of model constants 3,11 ). With a 
properly defined model error 3 ), the error in new regions of nuclei is about the same as in the 
region where the constants were adjusted. In the most recent investigations 11 ' 12 ) of the 1981 
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mass calculation in new regions of nuclei, the error for 351 new nuclei was only 6% larger than 
the error in the region where the model constants were adjusted. Furthermore, the error did not 
increase with distance from (3 stability. 

Also, many other nuclear-structure properties were successfully predicted by the model for 
nuclei far from stability 13-16 ). A special result of the 1981 mass calculation was the interpre- 
tation of certain spectroscopic results in terms of an intrinsic octupole deformation of nuclei in 
their ground state 1 ' 17 " 19 ). 

Here we present results of our new calculations of nuclear ground-state masses and defor- 
mations. Relative to the 1981 calculations we use an improved macroscopic-microscopic model, 
include additional shape degrees of freedom, extend the calculations to new regions of nuclei, 
and calculate a large number of additional nuclear ground-state properties. These additional 
properties will be published in a forthcoming article devoted to nuclear astrophysics 20 ). 

Specifically, we have improved the model in the following areas: 

• Our preferred macroscopic model is now the finite-range droplet model, which contains 
several essential improvements 4 ' 5 ) relative to the original droplet model 6 " 8 ). 

• The pairing calculations have been improved. Our pairing model is now the Lipkin-Nogami 
model 21-23 ). We also use an improved functional form of the effective-interaction pairing 
gap and an optimized pairing constant 9 ' 24 ). 

• An eighth-order Strutinsky shell correction is used. 

• The e zero-point energy is still added to the calculated potential energy to obtain the 
ground-state mass, but no 7 zero-point energy is added, since the method of calculation 
is not sufficiently accurate 11 ). 

• We minimize the ground-state energy with respect to €3 and e$ shape degrees of freedom, 
in addition to the 62 and £4 shape degrees of freedom considered previously. 

• Each ground-state shell-plus-pairing correction is based on single-particle levels calculated 
for the constants appropriate to the nucleus studied. Earlier, a single set of single-particle 
levels was used for an extended region of nuclei in conjunction with an interpolation scheme 
to improve accuracy. 

• The calculation has been extended from 4023 nuclei to 8979 nuclei, which now includes 
nuclei between the proton and neutron drip lines and superheavy nuclei up to A = 339. 

In the macroscopic-microscopic approach it is possible to calculate a large number of nuclear 
structure properties in addition to nuclear ground-state masses. These include the following: 

Even-multipole ground-state deformations: 



Quadrupole e deformation €2 

Hexadecapole e deformation 64 

Hexacontatetrapole e deformation e§ 

Related quadrupole j3 deformation /?2 

Related hexadecapole (3 deformation 04 

Related hexacontatetrapole (5 deformation (3q 
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Beta-decay properties: 

Q value of the (3 decay 
/3-decay half-life 

/3-delayed one-neutron emission probability 
/^-delayed two-neutron emission probability 
/3-delayed three-neutron emission probability 



Lipkin-Nogami pairing quantities: 

Neutron pairing gap A n 

Proton pairing gap A p 

Neutron number-fluctuation constant A2 n 

Proton number-fluctuation constant A2 P 

Odd-particle spins: 

Projection of the odd- neutron angular momentum along the symmetry axis $7 n 

Projection of the odd-proton angular momentum along the symmetry axis $7 P 



Alpha-decay properties: 

Q value of the a decay 
a-decay half-life 



Octupole properties: 

Octupole e deformation £3 

Related octupole (3 deformation /?3 

Decrease in mass due to octupole deformation AE3 



FRDM mass-related quantities: 

Spherical macroscopic energy 
Shell correction 
Pairing correction 
Microscopic correction 
Calculated mass excess 
Experimental mass excess 
Experimental uncertainty 
Discrepancy 

Calculated binding energy 

FRLDM mass-related quantities: 

Finite-range liquid-drop model microscopic correction 
Finite-range liquid-drop model mass excess 



1 l/2 
Pin 

P3n 



Qa 
1 l/2 
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Neutron and proton separation energies: 



One-neutron separation energy 5i n 

Two-neutron separation energy ,f>2n 

Three-neutron separation energy S$ n 

One-proton separation energy S\ p 

Two-proton separation energy Si p 



As mentioned above, we present here the calculated ground-state masses and deformations. 
Some of the remaining quantities will be presented in a forthcoming publication 20 ). 

In the next section we specify the macroscopic-microscopic finite-range droplet model in some 
detail. We discuss in particular the constants of the model, paying special attention to how to 
count the number of constants of a model. We present a summary of all constants in the model, 
including both those constants that have been determined from a least-squares adjustment to 
ground-state masses and fission-barrier heights and those that have been determined from other 
considerations. After our model has been specified, we discuss how it has been applied to the 
current mass calculation. 



2 Models 

In the macroscopic-microscopic method the total potential energy, which is calculated as a 
function of shape, proton number Z, and neutron number N, is the sum of a macroscopic term 
and a microscopic term representing the shell-plus-pairing correction. Thus, the total nuclear 
potential energy can be written as 

E pot (Z, N, shape) = E mSLC (Z, N, shape) + E S+P (Z, N, shape) (1) 

We study here two alternative models for E mSLC , given by Eqs. (|40|) and (|62|) below. The shell- 
plus-pairing correction is given by Eqs. (75) and ( |7q ) below. 



It is practical to define an additional energy, the microscopic correction E m i c , which is dif- 
ferent from the shell-plus-pairing correction E s+p . For a specific deformation e a , the latter is 
determined solely from the single-particle level spectrum at this deformation by use of Strutin- 
sky's shell-correction method 25 ' 26 ) and a pairing model. In contrast, the microscopic correction 
is given by 

Exnic\£aJ — ^s+p(^a) ^mac(fa) -E'mac^sphere) (2) 

This definition has the desirable consequence that the potential energy E po t of a nucleus at a 
certain deformation, for example the ground-state deformation e gs , is simply 

-^pot(£gs) — -^mic(Cgs) -^mac(Csphere) (3) 

However, the reader should note that the term microscopic correction is sometimes used instead 
for shell-plus-pairing correction. When results are presented it is usually E m { c that is tabulated, 
because it represents all additional effects over and above the spherical macroscopic energy. 
In practical calculations it is E s+p that is calculated. To obtain the total energy a deformed 
macroscopic energy term is then added to E s+p . 

There exist several different models for both the macroscopic and microscopic terms. Most 
initial work following the advent of Strutinsky's shell correction method used the liquid-drop 
model ' 28 ) as the macroscopic model. 
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The preferred model in the current calculations has its origin in a 1981 nuclear mass model li2 ), 
which utilized the folded- Yukawa single-particle potential developed in 1972 29,30 ). The macro- 
scopic model used in the 1981 calculation was a finite-range liquid-drop model, which contained 
a modified surface-energy term to account for the finite range of the nuclear force. The modi- 
fied surface-energy term was given by the Yukawa-plus-exponential finite-range model 31 ). The 
macroscopic part in this formulation does not describe such features as nuclear compressibility 
and corresponding variations in the proton and neutron radii. 

The droplet model 6-8 ), an extension of the liquid-drop model 27,28 ) that includes higher- 
order terms in A -1 / 3 and (N — Z)/A, does describe such features. However, in its original 
formulation the droplet model was very inaccurate for nuclei far from stability and also failed 
catastrophically 31 ) to reproduce fission barriers of medium-mass nuclei. These deficiencies led 
Myers to suggest that the surface-energy terms of the droplet model also be generalized to 
account for the finite range of the nuclear force. Thus, the Yukawa-plus-exponential model for 
the surface tension was incorporated into the droplet model. During this work it also became 
apparent that the description of nuclear compressibility was unsatisfactory, since the squeezing 
of the central density of light nuclei was overpredicted. The deficiency was serious because it 
starts to become important already at about A = 120 and becomes even more pronounced for 
lighter nuclei. To account for compressibility effects for light nuclei and for other higher-order 
effects an empirical exponential term was added 4 ' 5 ). 

The additions of the finite-range surface-energy effects and exponential term to the droplet 
model 5 ) resulted in dramatic improvements in its predictive properties, as summarized in the dis- 
cussion of Table A in Ref. 4 ). Mass calculations based on both the FRLDM 3 ) and the FRDM 4 ) 
were presented in the 1988 review of mass models in Atomic Data and Nuclear Data Tables. 
These calculations also used an improved pairing model relative to that used in the 1981 work. 
In the 1988 results the error in the FRDM was 8% lower than that in the FRLDM. 

However, there were two major unresolved issues in the 1988 calculations. First, there existed 
some deficiencies in the pairing model and the values of the constants that were used. Second, 
£3 and e% shape degrees of freedom were still not included, so deviations between calculated 
and measured masses due to the omission of these shape degrees of freedom were still present. 
Extensive investigations of pairing models and their constants have now been completed and 
resulted in an improved formulation of the pairing model 9 ). We have now also minimized the 
potential energy with respect to e-$ and e§ shape degrees of freedom. An overview of the results 
has been given in a paper on Coulomb redistribution effects 10 ). The FRDM, which includes 
Coulomb redistribution effects, is now our preferred nuclear mass model. 

Although the FRDM is now our preferred model, we also present results for the FRLDM 
for comparative purposes and for use in studies that assume constant nuclear density. We 
therefore specify below both models. Because several of the model constants are determined by 
least-squares-minimization of the model error, we start by defining model error. 



2.1 Model error and adjustment procedure 



In many studies the model error has been defined as simply the root-mean-square (rms) devia- 
tion, which as usual is given by 



rms 



1 n 

-E(^x P 



(4) 



Here M^ h is the calculated mass for a particular value of the proton number Z and neutron 
number N, and Mg xp is the corresponding measured quantity. There are n such measurements 
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for different N and Z. The choice (Q) is a reasonable definition when all the errors <7 exp associated 
with the measurements are small compared to the model error. However, for large cr* the above 
definition is unsatisfactory, since both the theoretical and experimental errors contribute to the 
rms deviation. The definition (Q) will therefore always overestimate the intrinsic model error. 

When the experimental errors are large, it is necessary to use an approach that "decouples" 
the theoretical and experimental errors from one another. This can be accomplished by observing 
that the calculated masses are distributed around the true masses with a standard deviation oth- 
There exist powerful statistical methods for determining the intrinsic model error o"th- The model 
error obtained in this way contains no contributions from the experimental uncertainties <7g Xp . 
To introduce such an error concept a new set of equations for determining model parameters and 
error were derived 3 ) by use of statistical arguments and the maximum-likelihood (ML) method. 
Here we generalize from the original assumption 3 ) e! h € N(0,o"th) that the theoretical error term 
e\ h is normally distributed with zero mean deviation from the true mass to e\ h S N(/ith,o"th) 
to allow for an error with a mean jj,^ that is different from zero and a standard deviation 
around this mean 12 ). This leads to the generalized equations 



" [M^ - (M* h + Mth *)] dMl h 



=1 u exp 



' , + <7 t h 2 " dp v 



0, v = l,2,...,m (5) 



^ [mi xp - (M{ h + m *)r - «, + ^h 2 *) _ 

2^ , ■ 2 o*\2 ~~ ' 

<=1 (^exp +°th Z ) 

^ [M- xp -(M- h + / x th *)] _ n 

i=l (^cxp +°"th^ ) 

where p u are the unknown parameters of the model. The notation o"th 2 * means that by solving 
Eqs. (|) and (@) we obtain the estimate a t h 2 * of the true <r t h 2 - Equation (||) is equivalent to 
minimizing S with respect to p v , where 

s _^ [ML p -(Mt h + y th *)f 

^ rri 2 4- n , 2 * 

i=l a exp + a th 

Thus, we are led to two additional equations relative to the usual least-squares equations 
that arise when model parameters are estimated by adjustments to experimental data under the 
assumption of a perfect theory, <7th = 0. For the FRLDM the least-squares equations (||) are 
linear, whereas for the FRDM they are non- linear. 

When the model contains a term cqA° that is strictly constant, Eq. (0) is identical to the 
member in Eq. (|5|) that corresponds to the derivative with respect to this constant. Thus, one 
should in this case put /i t h* = and solve only the remaining m+1 equations. One may therefore 
in this case characterize the error of the model in the region where the parameters were adjusted 
solely by the quantity a t h- In other cases one should solve the full set of equations. If ^ t h* is 
significantly different from zero the theory needs modification. Even if /x t h = in the original 
data region, it is entirely possible (although undesirable) that one obtains a mean error ^ t h* 
that is substantially different from zero when one analyzes model results for new data points to 
which the parameters were not adjusted. In this case the most complete characterization of the 
theoretical error requires both its mean /z t h and its standard deviation <7 t h around this mean. 

To allow for a single error measure that is similar to an rms deviation between the data and 
model we later also calculate the square root of the second central moment of the error term, 
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°th;//=0) m our studies of model behavior in new regions of nuclei. This quantity is obtained by 
setting //th* = when solving Eq. (^). In contrast to the rms measure, it has the advantage 
that it has no contributions from the experimental errors. 

Equations @-(@) constitute a system of m + 2 equations that are to be solved together. It 
is instructive to rewrite Eqs. (||) and ([?]) as 



1 



o~th 



El 



w 



17 E 



w. 



(M, 



exp 



fJ-th 



o~, 



exp 



(9) 



Ei 



E 



exp 



th/ 



where 



W; 



(^xp 2 + ^th 2 *) fe 



(10) 



(11) 



K = 2 (12) 

k„ = 1 (13) 

The unknowns // t h* and a t h 2 * can easily be determined from Eqs. (||) and ([To|) by an iterative 
procedure whose convergence is extremely rapid, requiring only about four iterations. An in- 
terpretation, not a proof, of Eq. @ is that the experimental error is "subtracted out" from the 
difference between the experimental and calculated masses. 

A common misconception is that one has to "throw away" data points that have errors that 
are equal to or larger than the error of the model whose parameters are determined. When a 
proper statistical approach, such as the one above, is used, this is no longer necessary. 

We will see below that the discrepancy between our mass calculations and measured masses 
systematically increases as the size of the nuclear system decreases. It is therefore of interest 
to consider that the mass-model error is a function of mass number A. A simple function to 
investigate is 

*th = ^7 (14) 



where c and a are two parameters to be determined. Whereas under the assumption of a 
constant model error one determines this single error constant from Eq. (S), we find that the 
ML method for the error assumption in Eq. (|l4|), with two unknowns, and assuming /i t h = 0, 
yields the equations 



E 



(MLp - M i 



th) 



a, 



exp 



+ 



At 



(7, 



exp 



+ 



At 



(15) 



n (ii4p-A4y 
E- 



'-'"exp 1 \ j^c 



i=l 



'-'exp 



A; 



(16) 



a*+l 



These equations are considerably more complicated to solve than Eq. (^). We solve them by 
minimizing the sum of the squares of the right members of Eqs. (^) and (|l6|). 
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2.2 Shape parameterizations 

The original parameterization of the folded- Yukawa single-particle model was the three-quadratic- 
surface parameterization 29 ' 32 ). It was designed to allow great flexibility in describing shapes late 
in the fission process. However, it is less suitable for describing ground-state shapes. 

To allow a better description of ground-state shapes and to allow close comparison with re- 
sults of Nilsson modified-oscillator calculations, we incorporated the Nilsson pertur bed-spheroid 
parameterization, or e parameterization, into the folded- Yukawa single-particle computer code 

in 1973 30,33,34) 

In our work here we use the e parameterization for all calculations related to ground-state 
properties. In our adjustment of macroscopic constants we also include 28 outer saddle-point 
heights of fission barriers. The shapes of these saddle points were obtained in a three-parameter 
calculation in the three-quadratic-surface parameterization in 197 3 33 ). 



2.2.1 Perturbed-spheroid parameterization 

The e parameterization was originally used by Nilsson 35 ) in the modified-oscillator single-particle 
potential. It was introduced to limit the dimensions of the matrices from which the single- 
particle energies and wave functions are obtained by diagonalization. This requirement leads to 
somewhat complex expressions for the nuclear shape. Here we employ its extension to higher- 
multipole distortions. For completeness we define it with axially asymmetric shapes 36-38 ) in- 
cluded, although this symmetry-breaking shape degree of freedom has not yet been implemented 
in the folded- Yukawa single-particle model. Note that the factor \J~^-\ is missing in front of the 
F 4 (7) function in Eq. (3) of Ref. 38 ). 

As the first step in defining the e parameterization a "stretched" representation is introduced. 
The stretched coordinates £, r], and ( are defined by 



c 



mujQ 



1 - ^e 2 cos (j + ^vr 



2 / 2 

1 - -e 2 COS (7 - -7T 



1 - 7^2 cos 7 



1/2 



1/2 



1/2 



(17) 



where Tiloo is the oscillator energy, e 2 the ellipsoidal deformation parameter, and 7 the non- 
axiality angle. It is then convenient to define a "stretched" radius vector p t by 



Pt = (e 2 + v 2 + c 2 ) 1/2 



(18) 



a stretched polar angle 9t by 



u = cos t>t 



Pt 



1 - -e 2 cos 7 



1 /1\ 1/2 
1 — g£2 cos 7(3 cos 2 6 — 1) + ( - J e 2 sin 7 sin 2 9 cos 2<f> 



1/2 



cos 9 (19) 
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and a stretched azimuthal angle 0t by 

2?/ 



v = cos 2</3 t 



1 + o £ 2 cos 7 



1X1/2 

cos 20 + ( - J e 2 sin 7 



(£ 2 + ti 2 ) 1 / 2 



1 /1\V2 
1 + g £2 cos 7 + I -1 e 2 sin 7 cos 20 



(20) 



In the folded- Yukawa model the single-particle potential is very different from that in the 
Nilsson modified-oscillator model. However, the definition of the e parameterization will be most 
clear if we follow the steps in the Nilsson model. The implementation in the folded- Yukawa model 
will then be simple. The Nilsson modified-oscillator potential is defined by 



V = hu; oPt 2 jl + 2eiPi(cos0 t ) 



2 1 / 8 ^ ^ 2 r 

- -e 2 cos7P 2 (cos6>t) + -e 2 sin 7 f — vr 



Y 2 {9 t ,<f> t ) + Y 



2 - 2 (0t,0 t ; 



+ 2e 3 P 3 (cos6l t ) + 2e 4 y 4 (cos/9 t ,cos20 t ) + 2e 5 P 5 (cos 6 t ) + 2e 6 P 6 (cos 6 t ) 



where 



kHloo 2l t ■ s + fi(l t 2 - < l t 2 >) 



Vi(u, v) = a 40 P 4 + \l^r a A2 {Y'i + Y^ 2 ) + a 44 (Y 4 4 + Y 4 



(21) 



(22) 



Here the hexadecapole potential V^u^v) is made dependent on 7 in such a way that axial 
symmetry is maintained when 7 = 0, 60°, —120°, and —60°, for mass-symmetric shapes and for 
€q = 0. This is accomplished by choosing the coefficients a 4 « so that they have the transformation 
properties of a hexadecapole tensor. However, this is achieved only for mass-symmetric shapes 
and for £q = 0. The e parameterization has not been generalized to a more general case. Thus 38 ) 



a 4 o 



-(5 cos 2 7 + 1) 

D 



a 42 



a 44 



12 



30 sin 27 



^\/70sin 2 7 



(23) 



It is customary to now assume that the shape of the nuclear surface is equal to the shape of 
an equipotential surface given by Eq. ( f2~l| ) . By neglecting the It • s and l t 2 terms and solving for 
p t and then using Eqs. (|l7|)-(|20|) to derive an expression for r in the non-stretched laboratory 
system we obtain 



KM) 



2 (2 
1 - -e 2 cos (7+ -vr 



2 (2 
1 - -e 2 cos I 7- -vr 



1 - -e 2 cos 7 



-1/2 



1 e 2 cos7 e 2 cos 7 + €2 I cos 7 -I — e 2 cos27 \u 



9 

l\ 1/2 
3 



e 2 sin 7 ^1 — -e 2 cos 7^ (1 — u 2 )v 



1/2 
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1 



1 - -e 2 cos 7^ (3u 2 - 1) + .. ( 



1\ 1/2 



€2 sin 7(1 — u 2 )v 
+ 2e 1 P 1 {u) + 2e 3 P 3 {u) + 2e 4 V 4 (u,v) + 2e 5 P 5 (u) + 2e 6 P 6 («) 



I-1/2 



(24) 



In the Nilsson model the starting point is to define the potential. After the potential has 
been generated the shape of the nuclear surface is deduced by the above argument. In the 
folded- Yukawa model the starting point is different. There, the equation for the nuclear surface, 
given by Eq. (24) in the case of the e parameterization, is specified in the initial step. Once 
the shape of the surface is known, the single-particle potential may be generated as described 
in later sections. 

The quantity ujo/ ujq is determined by requiring that the volume remain constant with de- 
formation, which gives 



1 

4-7T 



1 — 77^2 COS 7 + -7T 



1 — 77^2 COS 7 — -7T 



1 - -e 2 cos 7 



-1/2 









dO t a 


A sin O t 


1 


Jo Jo 







2 

— < 

3 



45 



1 -3/2 



+ 2exPx{u) + 2e 3 P 3 (u) + 2e i V i {u,v) + 2e 5 P 5 (u) + 2e 6 P 6 (n) 



(25) 



The above equation is derived by determining the volume inside the nuclear surface given by 
Eq. (p4|), with the integral / d?r inside the surface evaluated in terms of the "non-stretched" 
coordinates 9 and <f>. After a variable substitution one arrives at the expression in Eq. (p5|). 
The Legendre polynomials P\ occurring in the definitions of the e parameterization are defined 

by 



2 l l\ du l 

The first six Legendre polynomials are 

P (u) = 1 

Pi(u) 



1) 



1 = 0,1, 2, 



, 00 



(26) 



P2{U) 

Pa(u) 
P 5 (u) 
Psiu) 



II 



\(3u 2 - 1) 
i(5u 3 -3n) 



1 



(35u 4 - 30n 2 + 3) 



-(63n 5 - 70u 3 + 15m) 



1 



^(231m 6 - 315n 4 + 105m 2 - 5) 



The associated Legendre functions P ; m are defined by 



(1-M 2 ) m / 2 d l+m 
¥Tl du l + r 



1 = 0, 1, 2,... ,00; m = 0, 1, 2,...,l 



(27) 



(28) 
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The spherical harmonics are then determined from the relations 

■(2/ + 1) (/-m)P 1/2 



Yr(9,4>) = (-r 



PI" 1 (cos 6)e l 



4vr (l + m)\ 

Yr{eA) = (-) m Y l ~ m {6,4>) 

which yield for the functions used here 



m > 



(29) 
(30) 



Y 2 \9 

Y 2 \e 

Y£(9 

y 4 -\o 

Yi{6 
Y 4 ~ 2 (9 



15 
327 



sin 2 0e 2 ^ 



— sin 2 ee- 2 ^ 
32vr 



315 
512vr 

"315" 

512vr 

~45~ 

128vr 

"45~ 

1287T 



sin 4 Qe^ 



sin 4 6e-^ 



sin 2 0(7 cos 2 0- l)e 2 ^ 
sin 2 0(7 cos 2 0- l)e" 2 ^ 



(31) 



The sums 



SY 22 = Y 2 2 (0,0)+Kf 2 (M) 
SY 44 = Y 4 4 (9,<p)+Y^(e,4>) 

SY i2 = y 4 2 (0,0) + y 4 - 2 (0,^) 



(32) 



are required in the expression for the single-particle potential and in the corresponding equation 
for the nuclear surface. We obtain 



SY 22 
SY U 



sin 2 cos 2cj) 



(l-u 2 )v 



315 . 4 
sin 



128vr 



^ cos#= \/iE (i " u2)2(2 '' 2 " i) 



45 



SY 42 = \ — sin^(7cos 2 0-l)cos20 = J — (1 - u z )(7u z - l)v 



32vr 



2.2.2 Three-quadratic-surface parameterization 



45 



32vr 



(33) 



In the three-quadratic-surface parameterization the shape of the nuclear surface is specified in 
terms of three smoothly joined portions of quadratic surfaces of revolution. They are completely 
specified by 32 ) 



ai 2 - 


ai 2 
— {z 
ci 2 


-hf , 


h — c\ < z < z\ 




a 2 2 - 


a 2 i 
— {z 
c 2 z 


~hf , 


z 2 < z <l 2 + c 2 


(34) 


a 3 2 - 


G3 2 , 

— {z 

C3 2 


-hf , 


Z\ < z < z 2 
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Here the left-hand surface is denoted by the subscript 1, the right-hand one by 2, and the middle 
one by 3. Each surface is specified by the position l{ of its center, its transverse semiaxis ctj, and 
its semi-symmetry axis q. At the left and right intersections of the middle surface with the end 
surfaces the value of z is Z\ and Z2, respectively. 

There are nine numbers required to specify the expressions in Eq. ( |34| ) but three numbers 
are eliminated by the conditions of constancy of the volume and continuous first derivatives at 
z\ and z%. The introduction of an auxiliary unit of distance u through 



u 



(ai 2 + a 2 2 ) 



permits a natural definition of two sets of shape coordinates, 
coordinates <7j and three mass-asymmetric coordinates ctj by 

{h-h) 



(35) 

We define three mass-symmetric 



02 



Q3 
C3 2 



0"3 



a i 




a-2 



o 2 



a 3 = — y k (36) 

Ci z C 2 Z 

The coordinate a\ is not varied freely but is instead determined by the requirement that the 
center of mass be at the origin. 



a-2 



2.2.3 Conversions to (3 parameters 

A common parameterization, which we do not use here, is the (5 parameterization. However, 
since we want to present some of our results in terms of (3 shape parameters, we introduce the 
parameterization and a scheme to express shapes generated in other parameterizations in terms 
of (3 deformation parameters. In the (3 parameterization the radius vector r is defined by 

oo I 

r(9, 0) = R (l + E E Am^T) (37) 

1=1 m=—l 

where Rq is deformation dependent so as to conserve the volume inside the nuclear surface. When 
only axially symmetric shapes are considered the notation (3i is normally used for (3iq. Since the 
spherical harmonics Y™ are orthogonal, one may determine the (3 parameters corresponding to 
a specific shape in the e parameterization by use of 

/ r(e^)Yr(e,^dn 

(3lm = V^^j (38) 

J r(6,<f>)Y$(6,<l>)dn 
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where r is now the radius vector in the e parameterization, given by Eq. (|24|). This conversion 
equation is in fact valid for a radius vector r(8, cj)) defined by any parameterization. 

When the j3 parameters corresponding to a specific shape in the e parameterization are 
determined one should observe that higher-order (3 parameters may be non-zero even if higher- 
order e parameters are identically zero. For this reason, and because (3§ is not tabulated, the 
nuclear ground-state shape is not completely specified by the (3 parameters in the Table, whereas 
the shape is completely defined by the e parameters. 

2.3 Finite-range droplet model 

The finite-range droplet model, developed in 1984 5 ), combines the finite-range effects of the 
FRLDM 31,39 ' °) with the higher-order terms in the droplet model. In addition, the finite-range 
droplet model contains the new exponential term 

- CAe~^ Al/3 e (39) 

where C and 7 specify the strength and range, respectively, of this contribution to the energy 
and the quantity e is a dilatation variable given by Eq. ( ]49| ) below. The exponential term leads 
to an improved description of compressibility effects and is crucial to the substantially improved 
results obtained in the finite-range droplet model relative to the original droplet model. The 
necessity for this empirical exponential term, which is discussed extensively in Refs. 5,41 ), is 
clearly demonstrated in Refs. 5,41 ) and by the results obtained in Sec. 3.2 below. 

Most of our results here are based on the finite-range droplet model for the macroscopic 
term. Relative to the formulation given in Ref. 5 ), which unfortunately has numerous misprints, 
we use a new model for the average neutron and proton pairing gaps. The complete expression 
for the contribution to the atomic mass excess from the FRDM macroscopic energy is obtained 
after minimization with respect to variations in e and d, where 5 is the average bulk relative 
neutron excess given by Eq. ( |47|) below. One then obtains 
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E mac (Z,N, shape) = 

M n Z + M n N 



-ai + J5 2 - ^Ke 2 ) . 1 



+ a 3 A l / 3 B k 
+ a ^° 
Z 2 



mass excesses of Z hydrogen atoms and iV neutrons 

volume energy 

surface energy 



+ ci 



A 1 / 3 ' 
c 2 Z 2 A l ^B r 



ci 



A 1 / 3 



c 5 Z 



2 B W B S 



Z 2 



Ca (iV - Z) 



+ 1/1 + 



curvature energy 
energy 

Coulomb energy 
volume redistribution energy 

Coulomb exchange correction 
surface redistribution energy 

proton form-factor correction to the Coulomb energy 

charge-asymmetry energy 



1/A , Z and N odd and equal 
, otherwise 



Wigner energy 



+ < 



' + A p + A n - 5 np , Z and N odd 

+ A p , Z odd and N even 

+ A n , Z even and N odd 

k + , Z and N even 

a c iZ 2 - 39 



average pairing energy 



energy of bound electrons 



(40) 

This expression differs from the corresponding one used in our earlier calculations 5 ) only in 
the form of the average pairing energy appearing in the next-to-last term. One should note 
that after minimization the exponential term is present only implicitly in Eq. (HO) through its 
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presence in Eq. ( |49| ) below. For the average neutron pairing gap A n , average proton pairing gap 
A p , and average neutron-proton interaction energy <5 np we now use 9 ' 24 ' 42 ) 

A - r ™ c B s (AV\ 



'"mac-Ss 

h 



A P = (42) 



^ np = B~^ (43) 
The zero reference point for the pairing energy now corresponds to even-even nuclei rather than 
to halfway between even-even and odd-odd nuclei. 
The quantities c\, C2, C4, and C5 are defined by 



it. 

5 r 



ci 



C2 = 336 \J + K ]C1 



1/1 18 



2 



5 / 3 \ 2/3 

C5 = 64^ <*> 

In Eq. ([iOD we have kept only the first term in the expression for the proton form-factor 
correction to the Coulomb energy, so that /o is given by 

l/145\r p 2 e 2 

/0 = "8 I 48" J ~W (45) 

The bulk nuclear asymmetry 5 is defined in terms of the neutron density p n and proton density 
Pp by 

5 = (46) 

Pbulk 

and the average bulk nuclear asymmetry is given by 

3 ci Z B V B S \ V 9 J 1 B S 2X 



1 + 16 Q A 2 /3 Bl J / ( V 1 + 4Q AV3 5i J (47) 
The relative deviation in the bulk of the density p from its nuclear matter value po is defined 

by 

e = -i^2 (48) 
3 po 

and the average relative deviation in the bulk of the density is given by 

e = (ce^ Al/3 ~2a 2 ^ + ht + c^bA /K (49) 
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The quantity B\ is the relative generalized surface or nuclear energy in a model that accounts 
for the effect of the finite range of the nuclear force. It is given by 



A- 2 / 3 
8^ 2 r 2 a 4 



e -|r-r'|/o 

|r — r'l/a 



d 3 r d 3 r' 



(50) 



where the integration is over the specified sharp-surface deformed generating shape of volume 
V. Since the volume of the generating shape is conserved during deformation we have 



(51) 



where Rq is the radius of the spherical shape. The relative Coulomb energy B% is given by 



B-, 



15 A~ 5 / 3 
32vr 2 r 5 



d 3 r d 3 r' 



V r 



1 



1 + 



1 r 



2 Ode 



|r— r'l/ode 



(52) 



The quantities B\ and B% are evaluated for Rq = r^A 1 / 3 . However, in the finite-range droplet 
model the equilibrium value i?don of the equivalent-sharp-surface radius corresponding to the 
nuclear density is given by the expression 



R den = r A^ 3 (l + e) 



(53) 



Thus, the actual value of the nuclear radius is determined by the balance between Coulomb, 
compressibility, and surface-tension effects as expressed by Eq. (p9|). To calculate this balance it 
is necessary to know the response of the surface-energy and Coulomb-energy terms B\ and B3 
to size changes. To account for this response we introduce the quantities B2 and B4, which are 
related to the derivatives of B\ and B3 . These derivatives are evaluated numerically and during 
this evaluation the radius R of the generating shape is varied around the value r$A 1 / 3 . 

The quantity B2, which as mentioned above is related to the derivative of the relative gen- 
eralized surface energy B%, is defined by 



Bo 



1 










x 2 B^ 


2x~o 


_dx V 



(54) 



X=XQ 



with 



R , roA 1 / 3 

x = — and xq = 

a a 



(55) 



The quantity B4 is related to the derivative of the relative Coulomb energy B3 and is defined by 



£4 = -yo' 



' d 


fB 3 Y 


djj 


V y J. 



(56) 



y=yo 



with 



R , roA 1 / 3 
and y 



Ode 



(57) 
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For spherical shapes the quantities B\, B2, -B3, and B4 can be evaluated analytically. One 
obtains 

B{ 0) = l-* + (l + Xo )(2 + - + ^ 1 )e- 2x ° 



xo \ x x z 



Bf ] = 1 - (l + 2x - xo 2 ) e 



3 yo 2 



B<°> = 1 + 5 



15 21 3 A 9 7 , * _ 2 „ n 



8y 8y 3 4 V 2y Vo 2 2y 3 
3 15 63 3 / 2 12 32 42 21 » , . 

+ ~ 7—* + 7— + — 9+— 1 + — 1 + — ^1 e 

2/cT 2y 3 4y & 4 Vj/o 2/cT 2/0 J 2/0 2/0 5 



(58) 



The expression S3 for the relative Coulomb energy yields the energy for an arbitrarily shaped, 
homogeneously charged, diffuse-surface nucleus to all orders in the diffuseness constant Oden- 

The 

constants in front of the integrals for B\ and B3 are chosen so that B\ and B3 are 1 for a sphere 
in the limit in which the range constant a and the diffuseness constant ad en are zero, in analogy 
with the definition of the quantities B s and Bq in the standard liquid-drop and droplet models. 
The quantities B2 and B4, which are related to the derivatives of B\ and -B3, respectively, were 
introduced above to treat the response of the nucleus to a change in size, resulting from a finite 
compressibility. The shape-dependent quantities B s , B v , B w , B^, and B r , which are defined 7 ) 
in the standard droplet model, are given by 



B K 



B v = - 



A~ 2 / 3 
4-Kro 2 

15A- 4 / 3 
167r 2 ro 2 



/ dS surface energy 

Js 

J j ^- — — - — W^j d 3 rd 3 r' neutron skin energy 



225 A~ 2 fff 1 \ 2 

B w = — — ^ — ^ / / — W ) dS dS' surface redistribution energy 



647r 3 ro 6 J Js \|r — r'| / 

^-1/3 r / 1 l \ 

Bk = / I — — h — — ) dS curvature energy 

87rr Js \Ri R2J 

1575^4~ 7 / 3 f f ( 1 \ 2 

B x = — — — ^ — =— / / I -j — W ) d 3 r d 3 r' volume redistribution energy 



647r 3 ro 7 J Jv \ |r — r 



(59) 



W = / / ] — ~ — 77 d 3 r d 3 r' (60) 



where 

W = 

^itr® 6 J Jv | r 

is the average of the Coulomb potential and i?i and R2 are the principal radii of curvature. 
2.4 Values of FRDM macroscopic-model constants 

The constants appearing in the expression for the finite-range droplet macroscopic model fall 
into four categories. The first category, which represents fundamental constants, includes 1 ' 2 ) 
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Mh = 7.289034 MeV hydrogen- atom mass excess 

M n = 8.071431 MeV neutron mass excess 

e 2 = 1.4399764 MeV fm electronic charge squared 

One should note that for consistency we here use the same values for the fundamental constants 
as in our 1981 mass calculation 1 ' 2 ). Results of a more recent evaluation of the fundamental 
constants appear in Refs. 43 ' ). 

The second category, which represents constants that have been determined from consider- 
ations other than nuclear masses, includes 4 ~ 4 ) 



a e \ = 1.433 x 10 -5 MeV electronic-binding constant 

K = 240 MeV nuclear compressibility constant 

r p = 0.80 fm proton root-mean-square radius 

ro = 1.16 fm nuclear-radius constant 

a = 0.68 fm range of Yukawa-plus-exponential potential 

aden = 0.70 fm range of Yukawa function used to 

generate nuclear charge distribution 

The third category, representing those constants whose values were obtained from consider- 
ation of odd-even mass differences 9 ' 24 ' 42 ) and other mass-like quantities, are 



mac — 


4.80 MeV 


average pairing-gap constant 


h = 


6.6 MeV 


neutron-proton interaction constant 


W = 


30 MeV 


Wigner constant 


L = 


MeV 


density-symmetry constant 


«3 = 


MeV 


curvature-energy constant 



It should be noted that the final calculated mass excess is strictly independent of the value used 
for r mac . This constant affects only the division of the mass excess between a macroscopic part 
and the remaining microscopic correction. We will therefore not include r mac when we later 
count the number of constants in our mass model. It is the pairing constant r m ; c which enters 
the microscopic model that affects the mass excess. It will be discussed below. 

Since fi t h = in our case, Eqs. (||) and (||) can be solved with the experimental data set 
of 1654 masses with Z > 8 and N > 8 45 ) and 28 fission-barrier heights to determine the 
remaining macroscopic constants and the error of our model. Because it is now clear that the 
measurements of the masses of 31_34 Na that are listed in the 1989 midstream evaluation of 
Audi 45 ) are in error, we have made four revisions. For 31 ' 32 Na we use early results of mass 
measurements at TOFI 46 ). The final, slightly different values appear in Ref. 47 ). For 33 Na we 
use results of new measurements at GANIL 48 ). The data point for 34 Na is excluded. 

To present all the macroscopic model constants together we list them here but discuss their 
adjustment later. These constants are 
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a 2 
J 



Q 



16.247 
23.92 
32.73 
29.21 
0.0 
0.436 



60 
0.831 



MeV 
MeV 
MeV 
MeV 
MeV 
MeV 
MeV 



volume-energy constant 
surface-energy constant 
symmetry-energy constant 
effective surface-stiffness constant 
^4° constant 

charge-asymmetry constant 

pre-exponential compressibility-term constant 

exponential compressibility-term range constant 



7 



The pairing constant r m i c which enters the microscopic model is also determined in a least- 
squares minimization with the above 1654 masses, although no barrier heights were included in 
its determination. Once the value of r m j c had been determined the adjustment routines were run 
again, this time with barriers included, to yield the final values of the constants listed above. 
The value of r m - lc will be given in the section on microscopic constants. The resulting error in 
the FRDM is a th = 0.669 MeV. 

For completeness we also specify the mass-energy conversion factor used in the interim 1989 
mass evaluation. In this evaluation the relation between atomic mass units and energy is given 



Although a more recent value has been adopted 43 ' 44 ' 49 ), it is the above value, consistent with 
the 1989 interim mass evaluation 45 ), that should be used if our calculated mass excesses in MeV 
are converted to atomic mass units. 

2.5 Finite-range liquid-drop model 

In the present version of our model the contribution to the atomic mass excess from the FRLDM 
macroscopic energy is given by 




1 u = 931.5014 MeV 



(61) 
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E^{Z,N, shape) = 

MhZ + M n N mass excesses of Z hydrogen atoms and N neutrons 

— a v (\ — K v I 2 ^j A volume energy 

+ a s (l — K s I 2 ^j B\A 2 ^ surface energy 

+ a A° A energy 



Z 2 

+ c i — 77T^3 Coulomb energy 

A 1 ' 6 

z 4 / 3 

— C4 , Coulomb exchange correction 

A 1 / 6 

Z 2 

+ f(kff p )—r- proton form- factor correction to the Coulomb energy 

— c a (iV — Z) charge-asymmetry energy 

TI . I , T , I 1/i , Z and A/" odd and equal \ 
+ W \ + { ' ,. . Wigner energy 

\ 1 , otherwise / 

f + A p + A n - (5 np , Z and N odd 
+ A p , Z odd and TV" even 



+ < 



average pairing energy 



+ A n , Z even and TV" odd 

k + , Z and N even 

— a c \Z 2 - 39 energy of bound electrons 

(62) 

This expression differs from the corresponding one used in our earlier calculations 1,2 ) only in the 
form of the average pairing energy appearing in the next-to-last term. For the average neutron 
pairing gap A n , average proton pairing gap A p , and average neutron-proton interaction energy 
<5 np we now use 9 ' 24 ' 42 ) 

An = (63) 

A P = ^ (64) 

= ^273 ( 65 ) 

The zero reference point for the pairing energy now corresponds to even-even nuclei rather than 
to halfway between even-even and odd-odd nuclei. 
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In the above expressions the quantities c\ and C4 are defined in terms of the electronic charge 
e and the nuclear-radius constant ro by 



ci 



3e2 
5 r 

5/3 \ 2/3 
C4 = iUJ Cl 



(66) 



The quantity / appearing in the proton form-factor correction to the Coulomb energy is given 
by 

1 ,. -', 2 "il r, ■>■")- 1 -07 



1 r p 2 e 2 

where the Fermi wave number is 



f(k F r p ) 



145 
~48~ 



327 , x , 1527 
(fcprp) 2 + 



2880 

'QvrZX 1 / 3 1 



1209600 



(67) 
(68) 

The relative surface energy B s , which is the ratio of the surface area of the nucleus at the 
actual shape to the surface area of the nucleus at the spherical shape, is given by 

A ~2/3 



The relative neutron excess / is 



- ( 97tZ X 
N — Z 



N - Z 



47rro 2 



dS 



(70) 



The quantity B\ is the relative generalized surface or nuclear energy in a model that accounts 
for the effect of the finite range of the nuclear force. It is given by 



B, = 



A -2/3 



87r 2 r 2 a 4 

The relative Coulomb energy B3 is given by 



2 - 



|r-r'|/o 



d 3 rd 3 r' 



B3 



15 A-W 
32vr 2 r 5 



d 3 r d 3 r' 



V r 



1 + 



1 r 



2 <2dc 



-|r-r'|/o de 



For spherical shapes the quantities B\ and B3 can be evaluated analytically. With 



x 



A 1 / 3 



and y 



A 1 ' 3 



one obtains 



B 



(o) _ 



1 - 



B. 



(o) _ 



1 - 



5_ 

yo 2 



+ (l + x )(2 + — + 



x x < 



(71) 



(72) 



(73) 



1 



15 21 

+ 



%o 8yo 3 



9 



7 



7 



2yo Vo 2 2y ' ; 



-22/0 



(74) 



The expression B3 for the relative Coulomb energy yields the energy for an arbitrarily shaped, 
homogeneously charged, diffuse-surface nucleus to all orders in the diffuseness constant a^cn- The 
constants in front of the integrals for B\ and B% have been chosen so that B\ and -B3 are 1 for 
a sphere in the limit in which the range a and diffuseness a^en are zero, in analogy with the 
definition of the quantities B s and Bq in the standard liquid-drop model. 
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2.6 Values of FRLDM macroscopic-model constants 

The constants appearing in the expression for the finite-range liquid-drop macroscopic model 
fall into four categories. The first category, which represents fundamental constants, includes 1,2 ) 



Mh = 7.289034 MeV hydrogen-atom mass excess 

M n = 8.071431 MeV neutron mass excess 

e 2 = 1.4399764 MeV fm electronic charge squared 

The second category, which represents constants that have been determined from consider- 
ations other than nuclear masses, includes > 2 ) 



a e i = 1.433 x 10 -5 MeV electronic-binding constant 
r p = 0.80 fm proton root-mean-square radius 

ro = 1.16 fm nuclear-radius constant 

a = 0.68 fm range of Yukawa-plus-exponential potential 

Oden = 0.70 fm range of Yukawa function used to 

generate nuclear charge distribution 

The third category, representing those constants whose values were obtained from consider- 
ation of odd-even mass differences 9 ' 24 ' 42 ) and other mass-like quantities, are 



r mac = 4.80 MeV average pairing-gap constant 

h = 6.6 MeV neutron-proton interaction constant 
W = 30 MeV Wigner constant 

It should be noted that the final calculated mass excess is strictly independent of the value used 
for r mac . This constant affects only the division of the mass excess between the macroscopic 
part and the remaining microscopic correction. We therefore do not include r mac when we later 
count the number of constants in our mass model. It is the pairing constant r m ; c which enters 
the microscopic model that affects the mass excess. It will be discussed below. 

Since //th = in our case, Eqs. @) and (||) can be solved with the experimental data set of 
1654 masses with Z > 8 and N > 8 45 ) and 28 fission-barrier heights to determine the remain- 
ing macroscopic constants and the error of our model. To present all the macroscopic model 
constants together we list them here but discuss their adjustment later. These constants are 





= 16.00126 


MeV 


volume-energy constant 




1.92240 


MeV 


volume-asymmetry constant 


a s 


= 21.18466 


MeV 


surface-energy constant 


K s 


2.345 


MeV 


surface-asymmetry constant 


«0 


2.615 


MeV 


A constant 


Ca 


= 0.10289 


MeV 


charge-asymmetry constant 



The resulting error in the FRLDM is cr th = 0.779 MeV. 



2.7 Microscopic model 

The shell-plus-pairing correction E s+p (Z, N, shape) is the sum of the proton shell-plus-pairing 
correction and the neutron shell-plus-pairing correction, namely 

E S+P (Z, N, shape) = E^ p (Z, shape) + E™f (N, shape) (75) 
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We give here the equations for the neutron shell-plus-pairing correction. Completely analogous 
expressions hold for protons. We have 

Eg* (N, shape) = Effift (N, shape) + Eff (N, shape) (76) 

Both terms are evaluated from a set of calculated single-particle levels. As before, the shell 
correction is calculated by use of Strutinsky's method 25,26 ). Thus 

N 

E^ (N, shape) = ]T a - £ neut (N, shape) (77) 
i=l 

where ei are calculated single-particle energies and E neut (N, shape) is the smooth single-particle 
energy sum calculated in the Strutinsky method. The pairing correction is the difference between 
the pairing correlation energy and the average pairing correlation energy, namely 

Egg (N, shape) = E™^ (N, shape) - ££ e c ut (N, shape) (78) 

where BJ!^(N, shape) is given by Eq. ( |103| ) below and £p C c ut (iV, shape) is given by Eq. ( |110| ) 
below. For the pairing correction we now use the Lipkin-Nogami 21-23 ) version of the BCS 
method, which takes into account the lowest-order correction to the total energy of the system 
associated with particle-number fluctuation. 

The single-particle potential felt by a nucleon is given by 

V = Vy + V s . . + V c (79) 

The first term is the spin-independent nuclear part of the potential, which is calculated in terms 
of the folded- Yukawa potential 

fi(r) = -^-^ / I w , dV (80) 

47ra po t d Jv |r - r\/a pot 

where the integration is over the volume of the generating shape, whose volume is held fixed at 
|7ri? pot 3 as the shape is deformed. The potential radius i? po t is given by 

-Rpot = -Rden + ^den ~ E den / R den (81) 

with 

Rdcn = r A 1 / 3 (l + e) (82) 
The potential depth V p for protons and potential depth V n for neutrons are given by 

V p = V s + VJ (83) 

V n = V S - VJ (84) 

The average bulk nuclear asymmetry 5 appearing in Eqs. (|83| ) and ( p4|) and average relative 
deviation e in the bulk of the density appearing in Eq. (82) are given by the droplet model and 
thus depend on the values of the droplet-model constants. The FRDM macroscopic constants 
are determined in a nonlinear least-squares adjustment, which requires about 1000 steps to find 
the optimum constants. In principle, these constants should then be used in the determination 
of the single-particle potential, the potential-energy surfaces should be recalculated with the new 
constants, a new mass calculation should be performed, and a new set of macroscopic constants 
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should be determined, with this iteration repeated until convergence. Because the calculation of 
potential-energy surfaces is extremely time-consuming, only one iteration has been performed. 

Furthermore, in determining the single-particle potential we have used the following early 
forms 50 ) of the droplet model expressions for 5 and e: 

, ( T 3ci Z 2 \ . / 9 J 1 \ 

S =[ I+ 8Q^)/{ 1 + 4QA^) (85) 
- 1 202 + 1? + <**)/!<: (86) 



The range 

a pot = 0.8 fm (87) 

of the Yukawa function in Eq. ( |S0"| ) has been determined from an adjustment of calculated single- 
particle levels to experimental data in the rare-earth and actinide regions 34 ). It is kept constant 
for nuclei throughout the periodic system. 

The spin-orbit potential is given by the expression 



V s . . = -A 



h \ 2 <T • Wl X p 



,2m nuc c/ h 

where A is the spin-orbit interaction strength, m nuc is the nucleon mass, er is the Pauli spin 
matrix, and p is the nucleon momentum. 

The spin-orbit strength has been determined from adjustments to experimental levels in 
the rare-earth and actinide regions. It has been shown 1,14>34 ) that many nuclear properties 
throughout the periodic system are well reproduced with A given by a function linear in A 
through the values determined in these two regions. This gives 



for protons and 



A p = 6.0 y-^j + 28.0 = 0.025,4 + 28.0 = k p A + l p 



A n = 4.5 (— ) + 31.5 = 0.01875A + 31.5 = k n A + l n (90) 
V240/ 



for neutrons. 

Finally, the Coulomb potential for protons is given by 



V c (r)=e Pc (91) 

JV r — r 



where the charge density p c is given by 

* = (92> 

The number of basis functions used in our calculations is 

N has = 12 (93) 
The overall curvature of the basis functions is chosen to yield 

= C^/A 1 / 3 (94) 

with 



y cur 



(95) 
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2.8 Microscopic pairing models 

Because of its basic simplicity, the BCS pairing model 51-54 ) has been the pairing model of choice 
in most previous nuclear-structure calculations i' 2 - 29 - 55 ). However, a well-known deficiency of the 
BCS model is that for large spacings between the single-particle levels at the Fermi surface, no 
non-trivial solutions exist. In practical applications, these situations occur not only at magic 
numbers, but also, for example, for deformed actinide nuclei at neutron numbers N = 142 and 
152. By taking into account effects associated with particle-number fluctuations, the Lipkin- 
Nogami approximation 21-23 ) goes beyond the BCS approximation and avoids such collapses. 

In solving the pairing equations for neutrons or protons in either the BCS or Lipkin-Nogami 
model, we consider a constant pairing interaction G acting between N2 — N\ + 1 doubly degen- 
erate single-particle levels, which are occupied by 7V"i nt nucleons. This interaction interval starts 
at level N\, located below the Fermi surface, and ends at level N 2 , located above the Fermi 
surface. With the definitions we use here, the levels are numbered consecutively starting with 
number 1 for the level at the bottom of the well. Thus, for even particle numbers, the last 
occupied levels in the neutron and proton wells are N/2 and Z/2, respectively. 

The level pairs included in the pairing calculation are often chosen symmetrically around 
the Fermi surface. However, for spherical nuclei it is more reasonable to require that degenerate 
spherical states have equal occupation probability. This condition cannot generally be satisfied 
simultaneously with a symmetric choice of levels in the interaction region. We therefore derive 
the pairing equations below for the more general case of arbitrary N± and 7V 2 . 

In the Lipkin-Nogami pairing model 21-23 ) the pairing gap A, Fermi energy A, number- 
fluctuation constant A2, occupation probabilities v k 2 , and shifted single-particle energies e k are 
determined from the 2(N 2 — Aq) + 5 coupled nonlinear equations 



N 2 



N tot = 2j2v k 2 + 2(Aq - 1) 



(96) 



Vk 



2 = ^ 1 



k = N 1 ,N 1 + l,, 



Mo 



£fc = e k + (4A 2 - G)v k 2 , k = Ni, Ni + 1, . . . , N 2 



(97) 

(98) 
(99) 



where 



A 2 = « 



K k=N x J \k=N 1 J fc=JVi 




E2 2 
Uk V k 



N 2 

E 

k=m 



E4 4 
Uk V k 



(100) 



u k 2 = l-v k 2 , k = N 1 ,N 1 + l,...,N 2 (101) 
The quasi-particle energies E k of the odd nucleon in an odd-yl nucleus are now given by 22 ) 



Ek 



n oil/ 2 

(e fc -A) 2 + A 2 +A 2 , k = N 1 ,N 1 + l,...,N 2 



(102) 



P. Moller, J. R. Nix, W. D. Myers, and W. J. Swiatecki/ 'Nuclear Masses 26 

In the Lipkin-Nogami model it is the sum A + A2 that is identified with odd-even mass differ- 
ences 22 ). We denote this sum by Aln- 

The pairing-correlation energy plus quasi-particle energy in the Lipkin-Nogami model is 
given by 

N 2 ^2 q N 2 N 2 

Ep.c. = ( 2ufc2 ~~ nfc ) efc — ~Fi 2 ( 2ufc4 ~~ nfc ) ~~ 4 ^ 2 51 u k 2 Vk 2 + Ei6 odd ,N tot (103) 

k=N! fc=7Vi k=Ni 

where are the single-particle energies and n^, with values 2, 1, or 0, specify the sharp distri- 
bution of particles in the absence of pairing. The quasi-particle energy for the odd particle 
occupying level i is given by Eq. (102), and 6*^,^0* is unity if N toi is odd and zero if iVtot is 
even. 



2.9 Effective-interaction pairing-gap models 

In microscopic pairing calculations the pairing strength G for neutrons and protons can be 
obtained from effective-interaction pairing gaps Ac n and A Gp given by 9 ) 

= ^ (104) 

^ g p = ^ (105) 

The dependence of the pairing strength G on the corresponding effective-interaction pairing gap 
Aq is obtained from the microscopic equations by assuming a constant level density for the 
average nucleus in the vicinity of the Fermi surface. This allows the sums in the equations to 
be replaced by integrals. The average level density of doubly degenerate levels is taken to be 

p = \gW (106) 

where g is the smooth level density that is obtained in Strutinsky's shell-correction method 
and A is the Fermi energy of the smoothed single-particle energy 29,56 ). Thus, we can make the 
substitution 

N 2 



E f( e >< - A) =}► p / f(x)dx (107) 



where 

yi 



|iVtat + Ni - 1 
P 



-\N tot + N 2 

2/2 = — (108) 

P 



The gap equation (^7|) can now be evaluated for an average nucleus, with the result 
1 1„ f V2 dx 

1 



G 2 J y^ Jx 2 + A G 2 



2 P 



In Jy 2 2 + A G 2 + y 2 - In J Vl 2 + A G 2 + y x 



(109) 
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From this expression, the pairing strength G in the BCS model can be determined in any region 
of the nuclear chart. 

The same expression may also be used in the Lipkin-Nogami case, but some reinterpreta- 
tions are necessary. It is now the energies occurring in Eq. (|97| ) that are assumed to be 
equally spaced. These are not precisely the single-particle energies but are related to them 



by Eq. (|99D . Thus, in order for to be equally spaced, the single-particle energies must 



be shifted downward by the amounts (4A2 — G)vk . Since the occupation probability is 
approximately unity far below the Fermi surface and zero far above, the corresponding single- 
particle energy distribution is approximately uniform far above and far below the Fermi surface 
but spread apart by the additional amount 4A2 — G close to the Fermi surface. Although this 
decrease in level density near the Fermi surface is accidental, it is in approximate accord with the 
ground-state structure of real nuclei, since the increased stability associated with ground-state 
configurations is due to low level densities near the Fermi surface 24,56 ). 

In the Lipkin-Nogami model, it is the quantity A + A2 that is associated with odd-even 
mass differences, whereas in the BCS model it is A only that should be directly compared to 
the experimental data. This leads to the expectation that there is a related difference between 
A G N and A G CS , the effective-interaction pairing gaps associated with the LN and BCS models, 
respectively. Since we determine the constants of the model for A G N directly from least-squares 
minimization, it is not necessary to specify exactly such a relationship. However, the above 
observation is of value as a rough rule of thumb, and to remind us to expect that the effective- 
interaction pairing gaps in the BCS and LN models are of somewhat different magnitude. 

The expression for the average pairing correlation energy plus quasi-particle energy -E p . c . 
in the Lipkin-Nogami model is obtained in a similar manner as the expression for the pairing 
matrix element G. For the average pairing correlation energy plus quasi-particle energy in the 
Lipkin-Nogami model we then obtain 



p.c. 



1_ 

2 P 



(l/2 — G)[V2— \jvi 2 + A G 2 ) + ( w - G) ( yi + Jyi 2 + A G 2 



+ -(G-AX 2 )pA G 



tan 



if 2/2 



tan 



if Vi 



AO, 



odd,AT t ot 



A G J \Aq 

The expression for A2 for an average nucleus is fairly lengthy. It is given by 



where 



A 2 



G ( A — C 



4 \B-C 



(110) 



(111) 



.4 



mm 




y 2 2 + A G A 

yi 2 + A G 2 



B 



2^ 



A G 2 p 
16 



tan 



tan 
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2.10 Shell correction 

The Strutinsky shell-correction method 25,26 ) requires two additional constants, the order p and 
the range 7. The shell correction should be insensitive to these quantities within a certain 
range of values. Their values can therefore be determined in principle by requiring the plateau 
condition to be fulfilled. We have found that for heavy nuclei this condition is indeed fulfilled, 
with the shell correction for nuclear ground-state shapes insensitive to the values of these two 
constants. However, for light nuclei this is no longer the case. Here the shell correction may vary 
by several MeV for a reasonable range of values of the range 7. Moreover, the shell correction 
often does not exhibit any plateau. This probably indicates a gradual breakdown of the shell- 
correction method as one approaches the very lightest region of nuclei, where the number of 
single-particle levels is small. 

In the present calculation we choose 

p = 8 (113) 
for the order in the Strutinsky shell-correction method. The corresponding range 7 is given by 

7 = C SI Hu; B s (114) 

with 

C sr = 1.0 (115) 

and B s given by Eq. (|70|). This choice lowers the error of the mass model to 0.669 MeV from 
0.734 MeV obtained with the same range coefficient but no dependence on surface area in a 
sixth-order correction. 



2.11 Values of microscopic-model constants 

The constants appearing in the expressions occurring in the microscopic shell-plus-pairing cal- 
culation fall into four categories. The first category, which represents fundamental constants, 
includes 

?Tinuc = 938.90595 MeV nucleon mass 

he = 197.32891 MeV fm Planck's constant multiplied 

by the speed of light and divided by 2ir 
e 2 = 1.4399764 MeV fm electronic charge squared 

The electronic charge squared has already been counted among the macroscopic constants. 

The second category, which represents constants that have been determined from consider- 
ations other than nuclear masses, includes 1>2 ' 29 ) 



Ccur — 


41 MeV 


basis curvature constant 


= 


52.5 MeV 


symmetric potential depth constant 


Va = 


48.7 MeV 


asymmetric potential depth constant 


^den = 


0.82 fm 


potential radius correction constant 


Bden = 


0.56 fm 2 


potential radius curvature-correction constant 


Opot = 


0.8 fm 


potential diffuseness constant 


kp = 


0.025 


proton spin-orbit A coefficient 


l P = 


28.0 


proton spin-orbit constant 


k Q = 


0.01875 


neutron spin-orbit A coefficient 


In = 


31.5 


neutron spin-orbit constant 
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The third category, representing those constants whose values were obtained from consider- 
ation of mass- like quantities, are 

-^bas = 12 number of basis functions 

p = 8 order of Strutinsky shell correction 
C sr = 1.0 Strutinsky range coefficient 

The fourth category, representing those constants whose values were obtained from a least- 
squares adjustment simultaneously with the macroscopic constants of the FRDM, includes only 
one microscopic constant, namely 

r mic = 3.2 MeV LN effective-interaction pairing-gap constant 

In addition, the following droplet-model constants, which have been determined in an earlier 
study 50 ), are used in the expressions for the average bulk nuclear asymmetry 5 and average 



relative deviation e in the bulk density that are used to calculate V p , V n , and i?den in Eqs. (83) 
|), and fl82]), respectively: 

a\ = 15.677 MeV volume-energy constant 

a2 = 22.00 MeV surface-energy constant 

J = 35 MeV symmetry-energy constant 

L = 99 MeV density-symmetry constant 

Q = 25 MeV effective surface-stiffness constant 

K = 300 MeV compressibility constant 

ro = 1.16 fm nuclear-radius constant 



Use of these values in Eqs. (p5|) and (|86|) leads to 

T (N- Z)/A + 0.0112Z 2 /A 5 / 3 

6 = 1 + 3.15/AV3 < 116 ) 

O-l 47 noonT2 0.00248Z 2 

^=-^T7F + - 3305 + A ,/3 ( 117 ) 

One could in principle carry through the iterations discussed above to obtain a consistent set of 
droplet-model constants for the macroscopic part and for the single-particle potential, but the 
required computational effort would be extensive. However, the value of ro is precisely the same 
as that used in the macroscopic model. 



3 Enumeration of constants 

It is always of interest to have a clear picture of exactly what constants enter a model. Naturally, 
anyone who sets out to verify a calculation by others or uses a model for new applications needs 
a complete specification of the model, for which a full specification of the constants and their 
values is an essential part. Also, when different models are compared it is highly valuable to 
fully understand exactly what constants enter the models. Unfortunately, discussions of model 
constants are often incomplete, misleading, and/or erroneous. For example, in Table A of 
Ref. ) the number of parameters of the mass model of Spanier and Johansson 58 ) is listed as 
12. However, in the article 58 ) by Spanier and Johansson the authors themselves list in their 
Table A 30 parameters plus 5 magic numbers that are not calculated within the mass model 
and must therefore be considered parameters, for a total of at least 35 parameters. 
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Table 1: Constants in the FRDM. The third column gives the number of constants adjusted to nuclear 
masses or mass-like quantities such as odd-even mass differences or fission-barrier heights. The fourth 
column gives the number of constants determined from other considerations. 



Constants 


Comment 


Mass- like 


Other 


M H , M n , e 2 


Macroscopic fundamental constants 





3 


a, adon, K 


Macroscopic constants from considerations 
other than mass-like data 


n 
u 


a 
D 


L, a 3 , W , h 


Macroscopic constants obtained 

in prior adjustments to mass-like data 


4 


u 


aj., a2, J , Q, ao, 
C, 7, c a 


Macroscopic constants determined by 
current least-squares adjustments 


Q 
O 


U 


C - TX% c 


Microscopic fundamental constants 





2 


-^dcni ^den? Ccur; 
^p? ^p? ^n? *^pot 


Microscopic constants 





10 




Microscopic constants determined 

from considerations of mass-like quantities 


3 





* mic 


Microscopic constant determined by 
current least-squares adjustments 


1 





Oi, a 2 , J, K, L, Q 


Droplet-model constants that enter the single- 
particle potential (see discussion in text) 








Subtotals 




1G 


21 


Total 






37 



We specify here all the constants that enter our model, rather than just those that in the 
final step are adjusted to experimental data by a least-squares procedure. We also include such 
constants as the number of basis functions used and fundamental constants like the electronic 
charge and Planck's constant. 

3.1 Constants in the FRDM 

The discussion in the previous section allows us to enumerate the constants in the FRDM model 
in Table |l]. From this list we see that the macroscopic-microscopic method requires relatively 
few constants. One feature of the model gives rise to a small complication when counting the 
number of constants. Droplet-model constants occur also in the determination of the single- 
particle potential. However, a different set of constants is used here because, as discussed above, 
one does not know what the optimum values are until the calculation has been completed. In 
principle, the calculation should be repeated with the new droplet-model constants defining the 
single-particle potential until convergence is obtained. In Table |l] we have counted the number 
of constants as if this procedure had been carried out. 

However, since the droplet-model constants used in the present calculations are different 
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Table 2: Constants in the FRLDM. The third column gives the number of constants adjusted to 
nuclear masses or mass-like quantities such as odd-even mass differences or fission-barrier heights. The 
fourth column gives the number of constants determined from other considerations. 



Constants 


Comment 


Mass- like 


Other 


M H , M„, e 2 


Macroscopic fundamental constants 





3 


0>eh r Qi r pi 
a, ttden 


Macroscopic constants from considerations 
other than mass-like data 


u 


r. 



W , n 


Macroscopic constants obtained 

in prior adjustments to mass-like data 


o 

I 


u 


ao, c a 


Macroscopic constants determined by 
current least-squares adjustments 


6 





TlC, 7Tl nuc 


Microscopic fundamental constants 





2 


^s, ^a, -^-den; -^den, Ccur, 

^P; ^p? ^11, , ^pot 


Microscopic constants 





10 


-^*baS) Pi C sr , T m [ c 


Microscopic constants determined 

from considerations of mass-like quantities 


4 





ay, a 2 , J, K, L, Q 


Droplet-model constants that enter the single- 
particle potential (see discussion in text) 


3 





Subtotals 




15 


20 


Total 






35 



in the microscopic part and in the macroscopic part, different counting schemes could also be 
employed. Since the droplet-model constants used in the microscopic expressions are obtained 
from four primary constants 50 ) and nuclear masses were used only to give rough estimates of 
these constants, one may not wish to regard them as determined from mass-like quantities. One 
of the four primary constants is the nuclear radius constant ro, which has the same value as 
we use in our macroscopic model. Therefore, only three remain that could be considered as 
additional FRDM constants. With this classification scheme the number of constants adjusted 
to mass-like quantities remains 16 and the total number of constants in the model increases 
from 37 to 40. Alternatively, if we do count the three primary constants as adjusted to nuclear 
masses, the total number of FRDM constants is 40, while the number adjusted to mass-like 
quantities increases from 16 to 19. 

3.2 Constants in the FRLDM 

The constants in the FRLDM, which are either identical to or similar to the constants in the 
FRDM, are enumerated in Table |2[ We mentioned in the discussion of the FRDM constants 
that the six constants in the last line of Table [l] would converge to the values of the same 
constants listed earlier in the table after a sufficient number of iterations. In the FRDM these 
constants therefore need not be regarded as additional constants. In contrast, in the FRLDM 
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they must be regarded as constants obtained from adjustments to mass-like quantities. However, 
as mentioned in the discussion of the FRDM constants, these constants are all obtained from 
three primary constants, so we only include three in this category. 

4 Results 

4.1 Determination of ground-state shapes and masses 

The adjustment of constants in the macroscopic model is simplified enormously because the 
ground-state shape and fission saddle-point shape are approximately independent of the precise 
values of these constants when they are varied within a reasonable range 59 ). We therefore 
calculate the ground-state deformation with one set of constants and subsequently determine 
the various terms in the mass expression at this deformation. The constants of the macroscopic 
model can then be adjusted, with the nuclear shapes remaining fixed. 

A significant advantage of this approach is that the effect of new features can often be 
investigated without repeating the entire calculation from the beginning, which would take 
about 100 hours of CRAY-1 CPU time. For example, when we investigated different pairing 
models and determined the optimum value of the pairing constant, we needed to recalculate only 
the pairing-energy term for each of the 8979 nuclei in our study. Since we have in the initial part 
of the calculation determined ground-state shapes and stored the corresponding ground-state 
single-particle levels for all nuclei on disc, we need only read in the single-particle levels, do 
the pairing calculation, and readjust the model constants to obtain the effect of a new pairing 
model or new pairing-model constant. Such a study takes only about 20 minutes of CRAY-1 
CPU time. 

Our determination of mass-model constants and ground-state nuclear masses involves several 
steps. We first briefly list these steps and then continue with a more extensive discussion. 

1. Potential-energy surfaces are calculated versus e2 and €4. In this calculation, which was 
actually performed already in 1987, the FRLDM as defined in Ref. 3 ) is used, except that 
for the pairing calculations the BCS approximation is used instead of the LN approxima- 
tion. From these potential-energy surfaces the ground-state €2 and £4 deformations are 
determined. 

2. The ground-state energy is minimized with respect to €3 and also with respect to €q for 
fixed values of €2 and €4. 

3. When the resulting ground-state shapes have been determined, single-particle levels are 
calculated for each nucleus at the appropriate deformation and stored on disc. The shell- 
plus-pairing correction is also calculated and stored on disc at this time. The shell-plus- 
pairing correction is then available for use in the calculation of ground-state masses and 
in the determination of macroscopic-model constants. It is the only microscopic quantity 
required for the mass adjustment. 

4. Now that the ground-state shapes have been determined, the various shape-dependent 
functions that occur in the macroscopic energy are evaluated at each appropriate ground- 
state shape and stored on disc. 

5. Analogous steps to those above for masses are carried out also for 28 fission-barrier heights. 
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6. Least-squares adjustments are now performed, with the nuclear masses weighted 80% and 
the fission-barrier heights weighted 20%. The macroscopic-model constants are determined 
and the ground-state masses are calculated. 

7. Finally, when the ground-state shapes and masses and fission-barrier heights are known, 
other properties such as /3-decay half-lives, /3-delayed neutron-emission and fission proba- 
bilities, and Q values for a decay are calculated. 

For the major portion of the potential-energy-surface calculation we have chosen the following 
grid: 

e 2 = -0.50 (0.05) 0.50, e 4 = -0.16 (0.04) 0.16 (118) 

When the ground-state minimum is outside this grid we have used instead the expanded, but 
less-dense grid: 

e 2 = -1.0 (0.1) 1.0, e 4 = -0.28 (0.07) 0.28 (119) 

For large values of e 4 the nuclear shapes develop somewhat unnatural wiggles. These wiggles 
can be removed and the energy lowered by use of higher multipoles in the specification of the 
nuclear shape 30,60 ). We include in the first step of our calculations one higher multipole, namely 
€q. However, since in this step we want to consider only two independent shape coordinates, we 
determine e% at each value of 62 and e 4 by minimizing the macroscopic potential energy for 240 Pu. 
For heavy nuclei the value of e$ obtained in such a minimization is approximately independent 
of the nucleus considered. On the other hand, for very light nuclei minimization with respect 
to cq (and in some cases with respect to e 4 ) leads to values corresponding to unphysical shapes. 
These arise because if the distance across a wiggle on the nuclear surface is of the order of 
the range of the Yukawa-plus-exponential folding function, the nuclear energy increases very 
little but the Coulomb energy decreases strongly with increasing deformation. For €q we avoid 
this difficulty by minimizing the energy for 240 Pu, which is sufficiently large that also with €q 
distortions included the wiggles on the surface is larger than the range of the Yukawa-plus- 
exponential function. In the light region we avoid unphysical values of e 4 by including only a 
physical range of values in our grid. 

We use the single-particle states of the folded- Yukawa single-particle potential to calculate 
the shell-plus-pairing corrections at each grid point. Although the constants of the single-particle 
potential depend on Z and N, for the determination of the ground-state values of €2 and e 4 we 
use the same set of calculated levels for a region of neighboring nuclei, since it is too time- 
consuming to repeat the diagonalization for each value of Z and N. However, when the same 
levels are used for a moderately large region of nuclei, the shell correction for a magic nucleus 
calculated in this way may differ by 1 MeV or more from the shell correction calculated with 
the single-particle potential appropriate to that particular nucleus. 

To overcome this difficulty we proceed by first noting that most constants of the single- 
particle potential have been determined for nuclei close to line of /3-stability. Because of this 
and because the radius of the single-particle potential is one of its most important constants, we 
reduce the Z and N dependence of the constants of the microscopic model to an A dependence 
only. We next divide the nuclear chart into regions of suitable size, choosing for each region one 
set of single-particle constants. The regions are centered about the mass numbers A = 16, 20, 
40, 60, 80, 100, 120, 140, 160, 180, 200, 220, 240, 260, 280, 300, 320, and 340. The individual 
values of Z and N for each region are taken to be the closest integers corresponding to Green's 
approximation 61 ) to the line of /3-stability, namely 
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For each nucleus with a mass number different from one of these central mass numbers 
we calculate the microscopic corrections for two sets of constants. For example, nuclei with 
201 < A < 239 are included in the A = 220 calculation and nuclei with 221 < A < 259 are 
included in the A = 240 calculation. To determine E s+p for, say, a nucleus with A = 225 we 
linearly interpolate in terms of A between the result for A = 220 and the result for A = 240. 
We find that such an interpolation gives results that agree to within a few-hundred keV with 
those obtained with a single-particle potential appropriate to the specific nuclei concerned. 

Once the ground-state values of €2 and €4 are determined in this way, the shell-plus-pairing 
corrections are recalculated at these ground-state shapes with the exact single-particle poten- 
tial appropriate to each of the 8979 nuclei. The slight approximation made in calculating the 
potential-energy surfaces affects only the calculation of the shape and has a negligible effect on 
the final energy. 

After the ground-state €2 and £4 deformations are determined, we investigate the stability of 
the ground state with respect to €3 and e% shape degrees of freedom. Rather than simultaneously 
varying all shape degrees of freedom, we instead vary €3 and e& separately, with €2 and €4 held 
fixed at the values previously determined. When €3 is varied e$ is set equal to the value used 
in the original minimization with respect to 62 and 64, and when e§ is varied €3 is set equal 
to zero. The deeper of the two minima obtained in these two minimizations is selected as the 
ground state. The importance of the €3 and e§ shape degrees of freedom is discussed further 
in Ref. 10 ). Because surface wiggles should not become too small relative to the range in the 
Yukawa-plus-exponential function, the €q minimization is carried out only for nuclei with A > 60. 

After the ground-state shapes are determined, the shell-plus-pairing corrections and shape- 
dependent macroscopic functions are calculated and stored on disc. The programs that use this 
information to determine the macroscopic-model constants and calculate ground-state masses 
are then run. Although the least-squares adjustment is a nonlinear one, it takes only a few 
minutes to find the optimum constant set in an 8-constant variation and to calculate the final 
mass table. At this point it is relatively simple to investigate alternative model assumptions. 
As an example, we discuss the results of one such investigation concerning the effect of varying 
the microscopic pairing constant r m \ c . 

In our earlier pairing-model studies 9 ) we determined r m j c in Eqs. ( |104j ) and (105) by min- 
imizing the rms deviation between pairing gaps calculated in the LN model and experimental 
pairing gaps. An alternative possibility is to find r m j c by minimizing the error in the mass 
model. Because our change in the order of the Strutinsky shell correction does influence slightly 
the pairing calculations through the determination of G from the effective-interaction pairing 
gap Ac a small change in r m [c could in principle be required to obtain an optimum pairing 
calculation. In a study of how the model error a t h, the Lipkin-Nogami pairing gap Aln> the 
theoretical-mass pairing gap A t h.mass> and the neutron separation energy S n depend on r m j c , 
we first calculate the ground-state shell-plus-pairing corrections for several values of r m j c . For 
each of these values we then determine a set of macroscopic-model constants and generate a 
full-fledged mass table. In this process we also obtain microscopic pairing quantities and neu- 
tron separation energies and compare with experimental values. Recall that Aln is the sum 
of the pairing gap A and the number fluctuation constant A2 that occur in the Lipkin-Nogami 
equations. The pairing gap Ath.mass is determined from odd-even theoretical mass differences. 
The results are summarized in Table ||. 

Ideally, the minimum deviation should occur for all quantities at the same value of r m ; c , which 
is almost but not quite the case. As seen in Table g all minima are close to the mass-model 
minimum at r m ; c = 3.2 MeV. We therefore choose this value of r m i c for our microscopic pairing 
calculations. Experimental pairing gaps determined from odd-even mass differences contain large 
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Table 3: Determination of the pairing constant r, 







rms 


rms 


rms 


^*mic 


Mass model 


A LN 


Ath.mass 




(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


3.1 


0.6746 


0.1740 


0.2035 


0.409 


3.2 


0.6694 


0.1691 


0.2044 


0.411 


3.3 


0.6695 


0.1675 


0.2091 


0.416 


3.5 


0.6733 


0.1745 


0.2287 


0.432 



errors arising from non-smooth contributions to the mass surface other than pairing effects, for 
example, from shape transitions and gaps in the deformed single-particle level spectra. Since such 
contributions are equally present in A exp and A t h.mass, they should cancel out approximately 
in the difference between these two quantities, if the mass model errors were sufficiently small. 
Consequently, the other non-smooth contributions to the mass surface are not expected to affect 
an rms minimization of Ath.mass- It is therefore of interest to note that our chosen value of 
?~mic is intermediate between what would have been obtained from considering Ath.mass 

and Aln 

deviations. 

The FRDM, which includes Coulomb redistribution effects, is now our preferred nuclear 
mass model. Relative to the work described in Ref. 10 ), the following further improvements have 
been incorporated into the model. First, it was found that the 7 zero-point energy could not 
be calculated with sufficient accuracy in our current model. It is therefore no longer included, 
whereas the 62 zero-point energy is still retained. Second, we have also returned to our orig- 
inal choice of basis functions corresponding to 12 oscillator shells for all A values, instead of 
using somewhat fewer basis functions for lighter nuclei 10 ). Third, we now use an eighth-order 
Strutinsky shell-correction with range 7 = 1.0huJoB s instead of our earlier choice of a sixth-order 
Strutinsky shell correction with the same range coefficient but no dependence on surface area. 
The change in zero-point energy reduced the error in the calculated neutron separation energies 
from 0.551 MeV to 0.444 MeV and the error in the calculated masses from 0.778 MeV 10 ) to 
0.773 MeV. The second and third improvements further reduced the separation-energy error to 
0.411 MeV and the mass-model error to 0.669 MeV. The rms error for Ath.mass nas decreased 
in a similar manner as the error in S n . Although the effect of the mass- model improvements on 
Aln is small, the effect on A t h. m ass is dramatic. Relative to our earlier pairing calculation 9 ), 
the improvement is more than 20%. It is no accident that both S n and Ath.mass showed similar 
improvements. Both are determined from mass differences between nearby masses, and such 
differences dramatically improved when the inaccurate 7 zero-point energies were excluded from 
the calculations. The constants of the final model were presented in an earlier section. 

As seen in Table ^, the error in our mass model is now 0.669 MeV. We have also performed 
a mass calculation with the FRLDM as the macroscopic model and identical shell-plus-pairing 
corrections as in the FRDM calculation. For the FRLDM the corresponding error is 0.779 MeV, 
which is 16% higher. 

Figure Q shows the results of the FRDM calculation. As usual, the top part shows the dif- 
ferences between measured masses and the spherical macroscopic FRDM contributions plotted 
against the neutron number N, with isotopes of a particular element connected by a line. These 
experimental microscopic corrections are to be compared with the calculated microscopic cor- 
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rections plotted in the middle part of the figure. When the macroscopic and microscopic parts 
of the mass calculation are combined and subtracted from the measured masses, the deviations 
in the bottom part of the figure remain. The trends of the error in the heavy region suggest 
that this mass model should be quite reliable for nuclei beyond the current end of the periodic 



system. This has been made all the more plausible by simulations discussed in Sec. 4.3 on ex- 
trapability. When €3 and e§ shape degrees of freedom are included in the mass calculations, it 
becomes clear that the FRLDM, which does not treat Coulomb redistribution effects, is deficient 
in the heavy-element region, as is seen in Fig. |. Thus, our preferred mass model is now the 
FRDM, which includes compressibility effects and the associated Coulomb redistribution. 

4.2 Compressibility 

We have earlier 4 ) studied how the discrepancy between measured masses and calculated masses 
depends on the compressibility constant K and on the new exponential term. In this earlier 
investigation we used the 1984 version of the FRDM 5 ). We found that the minimum error 
occurred for K = 324 MeV. For this value the rms deviation between calculated and experi- 
mental masses was 0.666 MeV. For the conventional value K = 240 MeV the rms error was only 
marginally higher, namely 0.676 MeV. Because of the relative insensitivity of the rms error to 
the value of the compressibility constant K we retained 4 ) the historical value K = 240 MeV, as 
we also do here. 

It is of interest to investigate the sensitivity of the model error to K also in the current 
version of the model. In Fig. || the solid circles connected by a solid line show the theoretical 
error in the mass model as a function of 1/K. For each value of K the constants in the model 
are determined by minimizing the same weighted sum of barrier and mass errors that is referred 



to under point 6 in Sec. |4.l| . However, only the theoretical error in the mass model itself is 
plotted. The arrow at K = 243 MeV indicates the optimum value of K obtained when the 
compressibility coefficient is varied along with the other constants. Thus, in our current model 
we obtain in a least-squares minimization a compressibility coefficient that is close to the value 
K = 240 MeV that was adopted from other considerations, but the determination is clearly 
subject to a large uncertainty. 

We also investigate how the error for nuclei with N > 65 depends on K. This dependence 
is shown as solid squares connected by a long-dashed line. The model constants have the 
same values as obtained from the adjustments corresponding to the solid circles. Thus, no new 
adjustment is performed to this limited region of nuclei; we only investigate the behavior of the 
error associated with this region. 

Finally, we show as open circles connected by a short-dashed line the result obtained in an 
adjustment without the exponential term. Here the minimum of the weighted sum of mass and 
barrier errors occurs at K = 451 MeV. The minimum of the function plotted, which is the mass 
error only, occurs at a slightly higher value of K. 

The relatively low curvature of the solid curve shows that K cannot be reliably determined 
from an adjustment to nuclear masses. The conventional droplet model value K = 240 MeV is 
consistent with the result we obtain in a least-squares adjustment to masses and fission barrier 
heights, but from the adjustment alone one would not be able to rule out that K has some other 
value in the range from somewhat below 200 MeV to about 500 MeV. 

The long-dashed curve shows that heavy nuclei disfavor high values of the compressibility 
coefficient. This observation has been made earlier and was taken as evidence for a Coulomb 
redistribution effect 10 ). 

The short-dashed curve giving the results without an exponential term in the mass model is 
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moderately incompatible with a compressibility near 240 MeV and completely rules out a signif- 
icantly lower value. However, our preferred treatment of the compressibility is the formulation 
that includes the exponential term, in which treatment the restrictions on K are the much less 
severe ones given above. 



4.3 Extrapability 

One test of the reliability of a nuclear mass model is to compare deviations between measured 
and calculated masses in new regions of nuclei that were not considered when the constants of 
the model were determined to deviations in the original region. This type of analysis was used 
earlier by Haustein 62 ). However, we here considerably modify his approach. In addition to 
examining the raw differences between measured and calculated masses, we use these differences 
to determine the model mean discrepancy //th from the true masses and the model standard 
deviation <r t h around this mean, for new regions of nuclei. Whereas the raw differences do not 
show the true behavior of the theoretical error because errors in the measurements contribute 



to these differences, by use of the ideas developed in Sec. |24| we are able to estimate the true 
mean /x t h and standard deviation <7th of the theoretical error term e^h- 

Since our new mass model was developed only recently, we cannot test its reliability in new 
regions of nuclei because sufficiently many new data points are not available. Therefore, we 
have resorted to a simple simulation, in which we adjusted the constants in the model to the 
same experimental data set that was used in our 1981 mass calculation 1,2 ). Consequently, this 
calculation is not completely identical to the one on which Fig. [I] is based. The differences 
between the 351 new masses that are now measured 45 ) and the calculated masses are plotted 
versus the number of neutrons from ^-stability in Fig. |4j. We observe no systematic increase 
in the error with increasing neutrons from /3-stability. For the new region of nuclei the square 
root of the second central moment is 0.686 MeV, compared to 0.671 MeV in the region where 
the parameters were adjusted, representing an increase of only 2%. In contrast, mass models 
based on postulated shell-correction terms and on a correspondingly larger number of constants 
normally diverge outside the region where the constants were determined 11 ' 12 ). 

To study more quantitatively how the error depends upon the distance from /3-stability, we 
introduce bins in the error plots sufficiently wide to contain about 10-20 points and calculate 
the mean error and standard deviation about the mean for each of these bins by use of the 
methods described in Sec. |2.1| . The results for our 1981 finite-range liquid-drop model and for 
our 1992 finite-range droplet model, but adjusted only to the same data set as was used in our 
1981 calculation, are shown in Fig. For each model the central, light-gray band representing 
the original error region extends one standard deviation on each side of zero. The solid dots 
connected by a thick black line represent the mean of the error for nuclei that were not considered 
when the constants in the model were determined. The dark-gray area extends one standard 
deviation on each side of this line. The properties of the two models displayed in Figs. |5| are 
summarized in Table |4j. 

To test the reliability of the FRDM for extrapolation beyond the heaviest known elements we 
have performed a rather severe test in which we adjust the constants in the model only to data 
in the region Z, N > 28 and A < 208. There are 1110 known masses in this region compared 
to 1654 in the region Z, N > 8 used in our standard adjustment. Thus, about one third of all 
known masses are excluded, with nuclei removed from both ends of the region of adjustment. 
We then apply the model with these constants to the calculation of all known masses in our 
standard region and compare the results to our standard model in Fig. ||. The error for the 
known nuclei is now 0.745 MeV, compared to 0.669 MeV with our standard model adjusted to 
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Table 4: Comparison of errors of two different mass calculations. The errors are tabulated both for the 
region in which the constants were originally adjusted and for a set of new nuclei that were not taken into 
account in the determination of the constants of the mass models. The error ratio is the ratio between 
the numbers in columns 8 and 3. 



Original nuclei New nuclei 



Model rms cr t h JV lmc rms /x t h cr th cr t h ;M =o Error 

(MeV) (MeV) (MeV) (MeV) (MeV) (MeV) (MeV) ratio 



FRLDM (1981) 0.835 0.831 351 0.911 -0.321 0.826 0.884 1.06 

FRDM (1992) 0.673 0.671 351 0.735 -0.004 0.686 0.686 1.02 



all known nuclei. Although there is a noticeable increase of the error in the regions that were 
not included in the adjustment, an inspection of Fig. |6| indicates that the increased error in the 
heavy region is not due to a systematic divergence of the mean error, but rather to a somewhat 
larger scatter in the error. 

In our standard model the mass excesses of 272 110 and 288 110 are 133.82 MeV and 165.68 MeV, 
respectively. In our restricted adjustment we obtain 133.65 MeV and 166.79 MeV, respectively. 
Thus, although 288 110 is 80 units in A away from the last nucleus included in the restricted 
adjustment, the mass obtained in this numerical experiment is only about 1 MeV different from 
that obtained in the calculation whose constants were adjusted to nuclei up to 50 units in A 
closer to the superheavy region. Since our standard calculation is adjusted so much closer to 
the superheavy region than is the numerical experiment, we feel that it should be accurate to 
about an MeV in the superheavy region. Since models with and without Coulomb redistribu- 
tion energies often differ by considerably more, the masses of superheavy elements could provide 
very strong further confirmation of the existence of Coulomb redistribution effects. A suitable 
nucleus for such a test is 272 110. The FRDM, which includes Coulomb-redistribution effects, 
predicts a mass excess of 133.82 MeV for this nucleus, whereas the FRLDM, which does not 
include Coulomb-redistribution effects, predicts 136.61 MeV. 

Figure [l] shows that as the lighter region is approached the error gradually increases in a 
systematic way. We have explored this possibility by first determining the model error for limited 
regions of nuclei by use of Eq. (||). We select A = 25(25)250 as centerpoints of the regions and 
define each region to extend from .A cen t er — 24 to ^center + 25. The errors in these restricted 
regions are shown as solid circles in Fig. [7]. Since the trend of the error looks approximately like 
c/A a we have determined the parameters of this assumed error function by use of the maximum- 
likelihood equations tfHfy and (111). We find c = 8.62 MeV and a = 0.57. The error function 



corresponding to these parameters is plotted as a solid line. 
4.4 Fission barriers 

Calculated heights of the outer peak in the fission barrier are compared to measured values in 
Table [|. The results are also shown graphically in Fig. ||. Extensive fission studies based on 
earlier and current versions of the models discussed here are presented in Refs. 30 ' 31 ' 33 ' 63-68 ). 
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Table 5: Comparison of experimental and calculated fission-barrier heights for 28 nuclei. 



z 


N 


A 


Experimental 
barrier 
(MeV) 


Calculated 
barrier 
(MeV) 


Discrepancy 
(MeV) 


48 


61 


109 


34.00 


35.69 


-1.69 


66 


94 


160 


27.40 


27.88 


-0.48 


76 


110 


186 


23.40 


21.21 


2.19 




112 


188 


24.20 


21.07 


3.13 


80 


118 


198 


20.40 


19.16 


1.24 


84 


126 


210 


20.95 


21.81 


-0.86 




128 


212 


19.50 


19.69 


-0.19 


88 


140 


228 


8.10 


8.41 


-0.31 


90 


138 


228 


6.50 


7.43 


-0.93 




140 


230 


7.00 


7.57 


-0.57 




142 


232 


6.30 


7.63 


-1.33 




144 


234 


6.65 


7.44 


-0.79 


92 


140 


232 


5.40 


6.61 


-1.21 




142 


234 


5.80 


6.79 


-0.99 




144 


236 


5.75 


6.65 


-0.90 




146 


238 


5.90 


4.89 


1.01 




148 


240 


5.80 


5.59 


0.21 


94 


144 


238 


5.30 


4.85 


0.45 




146 


240 


5.50 


4.74 


0.76 




148 


242 


5.50 


5.25 


0.25 




150 


244 


5.30 


5.78 


-0.48 




152 


246 


5.30 


6.27 


-0.97 


96 


146 


242 


5.00 


4.24 


0.76 




148 


244 


5.00 


5.05 


-0.05 




150 


246 


4.70 


5.69 


-0.99 




152 


248 


5.00 


6.07 


-1.07 




154 


250 


4.40 


5.51 


-1.11 


98 


154 


252 


4.80 


5.31 


-0.51 



4.5 Ground-state masses and deformations 

In the Table we tabulate our calculated ground-state deformations in the e parameterization, 
the corresponding coefficients (3 in a spherical-harmonics expansion, the atomic mass excesses 
and microscopic energies calculated in both the FRDM and FRLDM, and experimental masses 
and associated errors that were used in the adjustment of model constants. 

To give an overview, most of our calculated quantities are plotted versus N and Z in the form 
of color contour diagrams. The calculated ground-state deformations 62, £3, £4, and e§ are shown 
in Figs. |9 12 , and the corresponding coefficients /?3, (3&, and are shown in Figs. 13 -fl6|- We 



observe some features that are by now well-known. For example, the quadrupole deformation 
parameter €2 increases by about 0.05 for each deformed region below the actinide region. Oblate 
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deformations occur in transition regions on the heavy side of most deformed regions. The 
hexadecapole deformation 64 is large and negative in the beginning of deformed regions and 
large and positive in the end of deformed regions. The coefficients fa, fa, and fa have the 
opposite sign from the corresponding e deformations, whereas fa has the same sign as €2 but is 
roughly 10% larger. 

The microscopic correction is plotted in Fig. 17. The familiar doubly magic regions around 
1 5oSn5o, 1 5oSn§2, and 2 g2Pt>i26 stand out clearly. The center of the superheavy region is located 
at 2 115i7g. The large negative microscopic correction originating in the superheavy region 
extends a significant distance towards the southwest and reaches into the deformed actinide 
region. It is these large, negative microscopic corrections that have made possible the extension 
of the known elements as far as fjjgMt]^. As is seen in Figs. 18 and [H], the largest effects of 63 
and 6q in experimentally accessible regions occur around 222 Ra and 252 Fm, respectively. 

In Fig. we show the discrepancy between experimental and calculated masses in the form 
of a contour diagram versus N and Z. Above N ~ 65 there are only a few nuclei with an error 
marginally larger than 1 MeV. The noticeable errors near Z = 40, N = 56 are probably related 
to the unique 14 ) shell structure in this region and the reinforcement of the N = 56 shell closure 
for proton number Z = 40 and proton numbers just below. Such proton-neutron interactions 
are not accurately described within any simple single-particle effective-interaction framework. 
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Figure captions 

Fig. |l] Comparison of experimental and calculated microscopic corrections for 1654 nu- 

clei, for a macroscopic model corresponding to the finite-range droplet model. 
The bottom part showing the difference between these two quantities is equiv- 
alent to the difference between measured and calculated ground-state masses. 
There are almost no systematic errors remaining for nuclei above N = 65, for 
which region the error is only 0.448 MeV. The results shown in this figure rep- 
resent our new mass model. 

Fig. |2] Analogous to Fig. [I], but for the finite-range liquid-drop model, which contains 

no Coulomb-redistribution terms. This leads to the systematic negative errors 
in the heavy region, which indicate that the calculated masses are systematically 
too high. 

Fig. m Relation between the compressibility coefficient K and the theoretical error in the 

mass model. Calculated values are indicated by symbols, which are connected 
by curves to guide the eye. In our standard FRDM, which is our preferred 
model, the theoretical error depends only relatively weakly on the compressibility 
coefficient in the range 200 MeV < K < 500 MeV, as is shown by the solid circles. 
Without the exponential term a relatively high compressibility coefficient would 
be required. The error in the heavy region, shown by the solid squares, indicates 
that especially heavy nuclei disfavor high values of the compressibility coefficient. 

Fig. H] Calculation to show FRDM reliability in new regions of nuclei. Here we use a 

smaller set of measured masses to determine the constants of the model than in 
the full calculation shown in Fig. ||. The errors for nuclei not included in the 
adjustment are displayed in this figure. The error is only 2% larger in the new 
region compared to that in the region where the constants were determined. The 
two largest deviations occur for 23 O and 24 O, which probably indicates that this 
region of light very neutron-rich nuclei is outside the range of model applicability. 
Proton number 8 is the lowest value of Z that we consider in this model. 

Fig. |5| Comparison of the error behavior for two models applied to new nuclei versus 

the number of neutrons from /3-stability. 

Fig. |6] Test of extrapability of the FRDM towards the superheavy region. The top 

part of the figure shows the error of the standard FRDM. In the lower part the 
error was obtained from a mass model whose constants were determined from 
adjustments to the restricted set of nuclei with Z, N > 28 and A < 208. In the 
light region of nuclei there is no noticeable divergence of the results obtained 
in the restricted adjustment. In the heavy region there is some increase in the 
spread of the error, but no systematic divergence of the mean error. Based on 
the more detailed discussion in the text we deduce that our calculated masses for 
the superheavy elements are accurate to about an MeV. 

Fig. [7| Error in the mass calculation as a function of mass number A. The theoretical 

error has been determined for limited regions throughout the periodic system. 
The error represented by each solid circle is based on nuclei in a region that 
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Fig. | 
Fig. | 



Fig. M 



Fie. 11 



Fig. \E 



Fie. P 



Fie. 14 



extends 24 mass units below the circle and 25 mass units above the circle. The 
points are well approximated by the function 8.62 MeV/A 0,57 . 

Comparison of experimental and calculated fission-barrier heights for 28 nuclei. 
Isotopes are connected by lines. 

Calculated ground-state values of |e2 1 for 7969 nuclei with N < 200. Oblate 
shapes are indicated with horizontal black lines. About 14 deformed regions 
stand out, bordered or partially bordered by blue lines corresponding to magic 
nucleon numbers. The magnitude of the deformation in the deformed regions 
increases by about 0.05 with successively lighter regions or as one goes from 
neutron-rich to proton-rich regions. Deformed regions above N = 184 usually 
have very low fission barriers, and should have spontaneous-fission half- lives that 
are too short to be detectable. 

Calculated ground-state values of €3 for 7969 nuclei with N < 200. Most nuclei 
in the investigated region inside the black line are stable with respect to mass- 
asymmetric octupole deformations; only 640 nuclei are unstable with respect to 
these deformations. The largest effects of experimental significance are centered 
around 2 fjRai34. Effects beyond N = 184 are of no experimental significance, 
since these nuclei are too short-lived to be observed. 

Calculated ground-state values of 64 for 7969 nuclei with iV < 200. Character- 
istically, the values are large and negative in the beginning of major deformed 
regions and large and positive in the end of major deformed regions. In accor- 
dance with this general trend, £4 is large and positive near the rock of stability 
in the vicinity of 272 110 near the end of the deformed "actinide" region. 

Calculated ground-state values of e§ for 7969 nuclei with N < 200. The behavior 
of €q is less regular than that of the lower, even multipole distortions. 

Calculated ground-state values of | /3a | for 7969 nuclei with N < 200, which have 
been obtained by use of the transformation ( |38| ) from the e deformations. Oblate 
shapes are indicated by horizontal black lines. Comments given in Fig. I also 
apply here. 

Calculated ground-state values of \fy\ for 7969 nuclei with N < 200, which 



have been obtained by use of the transformation (35) from the e deformations. 



Comments given in Fig. 10 also apply here. 



Fie. m 



Fig.H 



Calculated ground-state values of 04 for 7969 nuclei with iV < 200, which have 
been obtained by use of the transformation (|38| ) from the e deformations. Com- 



ments given in Fig. 11 also apply here, but note that the sign of f3± is opposite 
that of £4 when significant deformations develop. 

Calculated ground-state values of f3§ for 7969 nuclei with N < 200, which have 
been obtained by use of the transformation (|3"8|) from the e deformations. Com- 



ments given in Fig. |12| also apply here, but note that the sign of (3§ is opposite 
that of e% when significant deformations develop. 
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Fig. [17] Calculated ground-state microscopic corrections for 7969 nuclei with N < 200. 

Well-known doubly magic regions at M^jSnso, 1 |oSng2, and 2 g2Pt>i26 stand out 
clearly. The minimum in the superheavy region is offset somewhat from 298 114i84 
and is located instead at 294 115i7g. An interesting feature, also present in our 
first mass calculation 2,14 ), is the rock of stability at 2 QgMti63- 

Fig. [ll| Calculated ground-state octupole instability for 7969 nuclei with N < 200. Only 

640 nuclei exhibit any instability with respect to this shape degree of free- 
dom. The largest effect in the experimentally accessible region is —1.41 MeV 
for 22 jAci 33 . 

Fig. [n] Calculated ground-state hexacontatetrapole instability for 7969 nuclei with N < 

200. The instability is relative to the energy corresponding to the macroscopic 
equilibrium value of e§. The largest effect is —1.29 MeV for 2 ggEsi52. 

Fig. |2^ Discrepancy between measured and calculated masses. Above TV = 65 only a few 

discrepancies are marginally more that 1 MeV. There is a gradual increase of 
the error towards the light region. The large, fluctuating error near TV = 60 is 
probably due to deviations between our simple effective interaction and the true 
nuclear force. It is well-known that for Z ~ 40 there is a re-enforcement of the 
N = 56 sub-shell closure. Such effects cannot be described within the framework 
of a single-particle model. 
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EXPLANATION OF TABLE 

Table. Calculated Nuclear Ground-State Masses and Deformations, Compared to 
Experimental Masses Where Available 

Z Proton number. The mass table is ordered by increasing proton number. The 

corresponding chemical symbol of named elements is given in parenthesis. 

N Neutron number 

A Mass Number 

€2 Calculated ground-state quadrupole deformation in the Nilsson perturbed-spheroid 

parameterization 

€3 Calculated ground-state octupole deformation in the Nilsson perturbed-spheroid 

parameterization 

€4 Calculated ground-state hexadecapole deformation in the Nilsson perturbed- 

spheroid parameterization 

€q Calculated ground-state hexacontatetrapole deformation in the Nilsson perturbed- 

spheroid parameterization 

€g ym Calculated ground-state hexacontatetrapole deformation in the Nilsson perturbed- 

spheroid parameterization 

fa Calculated quadrupole deformation of the nuclear ground-state expressed in the 

spherical-harmonics expansion (|3^) 

/?3 Calculated octupole deformation of the nuclear ground-state expressed in the 

spherical-harmonics expansion ([}?]) 

/?4 Calculated hexadecapole deformation of the nuclear ground-state expressed in 

the spherical-harmonics expansion ( |37| ) 

{3q Calculated hexacontatetrapole deformation of the nuclear ground-state expressed 

in the spherical-harmonics expansion ( |37| ) 

-Emic Calculated ground-state microscopic energy, given by the difference between the 

calculated ground-state atomic mass excess and the spherical macroscopic energy 
calculated from Eq. (|40|), in our preferred model, the FRDM 

M t h Calculated ground-state atomic mass excess, in our preferred model, the FRDM 

M exp Experimental ground-state atomic mass excess in the 1989 midstream evaluation 

of Audi ), with 4 revisions 

cr ex p Experimental error associated with the ground-state atomic mass excess in the 

1989 midstream evaluation of Audi 45 ), with 4 revisions 

E^ c Calculated ground-state microscopic energy, given by the difference between the 

calculated ground-state atomic mass excess and the spherical macroscopic energy 
calculated from Eq. <M), in the FRLDM 



M t ^ L Calculated ground-state atomic mass excess, in the FRLDM 
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